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IOTROttlCTION 

Knowledge  of  the  factors  which  contribute  to  pressure  drop  in  packed  beds 
finds  many  applications.  Motion  of  ground  waters,  petroleum,  and  natural  gas 
through  rocks,  soil,  and  sands  recuire  a  knowledge  of  the  same  laws  that 
govern  flow  through  packed  beds,  seepage  under  dams,  the  permeability  of 
concretes,  and  surface  area  or  density  of  many  industrial  materials  are 
determined  by  application  of  these  laws.  Direct  application  to  chemical 
engineering  is  found  in  filtration,  distillation,  absorption,  fluidized  or 
packed  bed  catalytic  operations,  and  drying  of  solid  particles. 

Several  attempts  have  been  made  to  describe  packed  beds  adequately  so 
that  pressure  drop  could  be  predicted.  These  investigations  have  success- 
fully answered  problems  of  limited  scope,  but,  none  have  resulted  in  suffi- 
ciently general  conclusions  to  allow  the  extension  of  existing  information 
to  new  packing  materials  of  novel  geometrical  nature. 

This  investigation  was  initiated  in  order  to  isolate  a  series  of  packed 
bed  variables  that  could  be  used  as  a  criterion  for  predicting  pressure  drop 
in  any  randomly  packed  bed.  Special  attention  was  given  to  amplifying  the 
effect  of  variables  that  had  already  been  recognized  as  important,  isolation 
of  the  relevant  factors,  and  a  general  formulation  of  these  results  that 
would  lead  to  accurate  prediction  for  pressure  loss  through  any  packing 
material. 

Variables  Y.'hich  Have  Been  Considered 

The  factors  which  affect  pressure  drop  can  be  segregated  into  two 
groups.  The  fluid  and  empty  column  variables  constitute  one  group.  The 
packed  bed  variables  constitute  another. 


The  role  of  fluid  variables  has  long  been  understood  sufficiently  so 
that  they  need  not  be  analysed  in  the  discussion.  It  is  sufficient  to 
select  nomenclature  for  these. 

Table  1.  Nomenclature  for  fluid  and  empty  column  variables. 

g  -  gravitational  constant  converting  weight  units  to  force  units. 

L^  =  depth  of  the  packed  zone. 

AP  ■  pressure  loss  due  to  frictional  resistance  across  L^«  Units 
of  weight  divided  by  area. 

U0  3  velocity  based  upon  the  empty  column. 

/A     =*   absolute  viscosity  of  the  fluid. 

P  ■  mass  density  of  the  fluid. 


Table  2  includes  basic  items  which  have  been  given  prior  consideration. 
More  complex  variables  are  not  included  since  each  different  investigator 
has  grouped  the  items  to  suit  his  theory  or  needs.  Expression  of  these  com- 
plex terms  has  been  avoided  for  the  sake  of  simplicity  of  understanding  the 
subject  as  a  whole. 

The  measurements  used  to  describe  packed  columns  as  treated  in  the 
literature  are  given  in  Table  2. 


Table  2.  Nomenclature  for  packing  variables  needed  for  discussion  of 
literature . 


D_  *  nominal  diameter  of  the  packing  unit. 

Ds  «  diameter  of  a  sphere  having  the  same  volume  as  the  packing 
unit. 

D+  ■  column  diameter. 

e   ■»  height  of  an  element  of  surface  roughness. 

r   *  ratio  between  distance  traversed  by  a  fluid  and  the  length 
of  the  column. 

s   ■  surface  of  the  packing  unit. 

S   ■  total  surface  of  the  packing* 

sg  ■  surface  of  a  sphere  having  the  same  volume  as  the  packing 
unit. 

S+  ■  surface  area  of  the  column. 

V  *  total  free,  or  void,  volume  in  the  packed  zone. 

v  *  volume  of  the  packing  unit. 

Vp  ■  total  volume  occupied  by  the  packing. 

V^  »  V  +  Vp  ■  volume  of  the  empty  column. 


Some  Important  Concepts 

One  very  basic  concept  of  the  problem  of  fluid  flow  has  become  classic 
to  the  representation  of  pressure  drop  in  packed  beds.  Dimensional  analysis 
based  upon  the  assumption  that  pressure  drop  is  a  function  of  fluid  density, 
fluid  viscosity,  fluid  velocity,  hydraulic  radius  or  diameter,  and  surface 
roughness  leads  to  results  that  parallel  the  development  of  the  Fanning 


equation  which  is  concerned  with  flow  through  channels.  The  result  is 

dP  .  (P\P     ifvPV,   e"\  (1) 

dL   6"^5"  P[  ~    B) 

where 

dP  ■  gradient  of  frictional  pressure  loss  along  the  actual  path 
dL   of  flow, 

O  -  constant  factor  determined  from  the  geometrical  nature  of 
the  available  flow  path, 

U  ■  mean  fluid  velocity  along  the  actual  path  of  flow, 

E  ■  the  equivalent  hydraulic  diameter  or  radius  of  the  bed, 

other  terms  are  as  in  Tables  1  or  2, 

aarly  developments  did  not  recognize  the  influence  of  o  and  *.,  The 

v 

efforts  of  each  individual  to  express  these  four  variables  through  consider- 
ation of  those  in  Table  2  and  the  geometrical  nature  of  the  packing  unit  has  con- 
stituted all  theoretical  developments  concerned  with  flow  through  packed  beds. 
According  to  the  results  of  this  investigation  a  satisfactory  expression 
for  hydraulic  radius  was  perhaps  first  introduced  by  Blake  (1)  in  1922,  Blake 
did  not  include  this  development  in  his  paper,  but  his  results  were  identical 
to  those  of  Carman  (£)  who  chose  to  solve  for 

volume  of  duct 


m 


surface  contacted  by  the  fluid 
This  was  carried  even  closer  to  the  concept  of  hydraulic  radius  as  it  was 
applied  to  flow  through  channels  by 


■ 


*  ^oss-section  of  duct  .  cross-section  x  length  of  duct 
perimeter  of  duct  ptrimeter  x  length  of  duct 

volume  of  duct 


surface  of  duct 


s 


;  the  mode  of  development,  it  is  certainly  trae  that  this  quantity 
is  a  measure  of  the  distance  between  surfaces  that  contact  the  fluid.  It  has 
been  found  convenient  to  write 

m  -  V/Sp  (2) 

No  sei'ious  difficulty  would  be  encountered  by  including  column  surface  with 
packing  surface,  a  step  -which  appears  logical.  However,  a  certain  degree  of 
mathematical  simplicit,,  ie  gained  by  use  of  equation  (2), 

Prediction  of  pressure  drop  through  a  wide  variety  of  packing  materials 
and  packed  bed  variables  was  finally  accomplished  by  introducing  a  concept 
which  might  be  termed  hydraulic  width  of  the  packing  unit.  In  any  event,  the 
new  term  represents  the  width  of  the  barrier  that  must  be  circumvented  by 
the  fluid  just  as  ,Tm"  represents  the  width  of  the  path  available  for  flow. 
A   definition  was  framed  so  that  the  following  could  be  approximated  mathemati- 
cally, 

surface  contacting  the  fluid 

yf    m  i  i  i    i     ii  i     n 

boundry  of  obstructing  surfaces 
To  this  end  the  packing  perimeter,  C_,  was  defined  ass  the  locus  of  tangent 

points  to  the  packing  that  would  be  generated  by  a  line  which  moved  through- 
out the  packed  bed  remaining  oriented  parallel  to  the  column  wall.  It  was 
then  possible  to  include  this  mathematical  expression  for  w: 

*  -  Sj/Cp  (3) 

By  combining  (2)  and  (3)  a  very  useful  measure  of  the  distortion  of  flow 

path  was  obtained  and  retained  in  the  form 

w 

—  •  index  to  degree  of  distortion  of  flow  path  and  other  effects 

of  the  packing. 
It  is  not  clear  at  present  whether  w/m,  w(l  *  S+/Sp)/m,  or  w(l  +  0,6Sj/Sp)/m 

will  be  proved  most  useful.     The  term,   (1  ♦  0,6S|/Sp),  was  discovered  as  a 


very  excellent  factor  for  accommodating  the  confining  surface  of  a  circular 

column;  its  development  will  be  treated  later  since  it  is  of  secondary 

importance • 

Plate  I,  Figs.  1,  2,  and  3  represent  the  way  in  which  1   indexes  the 

m 

degree  to  which  the  fluid  stream  is  disturbed.     Figure  1  illustrates  a  bed 
of  spheres,   and  the  arrows  suggest  the  path  of  flow  through  the  bed.     Figure 
2  is  constructed  approximately  to  scale  so  that  a  bed  of  packing  units,  which 
is  composed  of  randomly  suspended  circular  plates,   hss  the  same  values  of 
"m"  and  "w"  as  does  the  bed  of  spheres.     According  to  the  results  of  the 
investigation,   the  pressure  drop  through  both  beds  will  be  identical  if 
U0, /°,   and/d   are  identical.     Figure  3  represents  a  bed,    similar  to  the  one 
in  Fig.  2,   that  has  the  same  value  for  "m"  but  fewer  units  of  greater  HwH 
are  contained.     It  is  predicted  that,  if  U0,  /°,   and^/  are  the  same  as 
before,   the  pressure  drop  through  this  bed  will  be  greater  than  that  for 
either  of  the  other  beds.     On  the  other  hand,  wire  packing  producing  the 
same  value  for  "m"  should  allow  very  low  pressure  drop. 

The  illustration  was  not  tested  experimentally,  but  the  principle  was 
repeatedly  tested  by  reference  to  pressure  drop  through  packed  beds  of 
widely  differing  properties  which  were  approximately  equivalent  to  those  in 
the  example. 


no         .1 

Fig.  1.  Flon  through  spheres  of  diameter  IX  producing  A?1  for  ra-i,  ^1,  /*  1, 

Fig,  2,  Flow  thro.    L*tea  cf  diameter  2   x  I  should  produce^-,  for  bh,  ^1, 

Fig.  3»  Flow  through  pistes  of  diameter  I  x  Di  should  produce  AP  >AP,  for  m\9 
^1,^1,  (U0)!. 


PLATE  I 


figure 


figure  2 


figure  3 


LI- 

The  need  for  an  understanding  of  the  factors  which  contribute  to  pressure 
drop  in  packed  beds  was  realized  as  early  as  l863«  The  earlier  theories 
utilised  the  assumptions  that  a  packed  bed  was  comprised  of  a  series  of  ducts 
that  possessed  a  total  sectional  area  equal  to  the  area  that  would  be  inter- 
sected by  a  plane  passing  through  the  bed,  and  that  the  surface  of  the  walls 
of  these  ducts  was  equivalent  to  the  surface  of  the  packing  material.  Inves- 
tigators who  have  considered  one  or  two  different  packing  materials  have  been 
relatively  successful.  Those  who  have  considered  a  large  number  of  packing 
.-.aterials  have  accepted  serious  discrepancies  or  have  resorted  to  extremely  com- 
plex and  unjustifiable  empiricisiTis, 

In  1922  Blake  (1)  analysed  flow  through  several  beds  of  glass  and  clay 
Raschig  rings.  He  applied  the  principles  of   eouation  (1)  with  a  relative 
degree  of  success.  In  developing  a  friction  factor  and  Reynold* s  number,  he 
assumed  that  velocity  should  be  in  excess  of  the  superficial  velocity,  U0, 
by  the  ratio  V^/V  and  that  hydraulic  radius  should  be  expressed  as  in  equation 
(2).  The  ratio,  V/V^,  had  been  previously  proven  identical  to  the  fraction 
of  column  cross-section  that  wss  not  intersected  by  the  packing  material. 
Blake  represented  the  friction  factor  by  APgVvItUo  SpVtf°and  the  Reynold's  num- 
ber by  /°U0Vt^aSp.  These  terms  were  plotted  on  logarithmic  coordinates.  He 
found  that  a  single  line  represented  results  of  tests  with  the  glass  packing 
but  that  lower  values  of  the  friction  factor  were  obtained  for  the  clay 
packing. 
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Flow  through  beds  comprised  of  lead  shot  was  studied  by  Burke  and 
Plummer  (U)«  They  utilized  variables  equivalent  to  those  of  Blake  and  also 
represented  their  results  on  logarithmic  coordinates.  These  tests  illus- 
trated the  manner  in  which  fa   of  equation  (1)  depends  upon  the  Reynold* s  num- 
ber representation.  The  groupings  that  were  used  were  equivalent  to  those 
used  by  Blake.  Figure  3>  illustrates  these  results.  The  region  where 
AP/JL^  is  proportional  to  U0  is  laminar  flow  region.  In  the  turbulent  region 
AP/L-fc  is  proportional  to  U0  •  The  intermediate  zone  of  transition  from 
laminar  to  turbulent  flow  is  peculiar  to  packed  beds.  A  single  line  repre- 
sented the  results  with  a  satisfactory  degree  of  accuracy. 


friction 
factor 


Reynold's      number 


Fig.  5.     Results  of  Burke  and  Plummer  (h) • 
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Literature  up  to  shout  ]      8  thor        rvoyed  by     m   (5).  T'his 
survey  included  much  information  that  did  not  reach  journals  in  English  print. 
He  found  that  laminar  flow  could  be  represented  conveniently  by  the  equation 

It     ^ 
The  proportionality  factor  included  considerations  of  o  in  equation  (1)  from 

the  standpoint  of  the  nature  of  a  sphere  and  an  emperical  -wall  factor  correc- 

tion,  (1  ♦  St/Sp),  which  yielded  AP/l^  »  6  (i  +  ^/Bp)U0Sp2Vt/KV3.  Sam 

results  of  tests  by  Pirie,  given  to  Carman  in  a  private  communication,  supply 

information  for  beas  packed  with  cubes  and  prisms.  Pirie  was  quoted  to  have 

worked  "entirely  in  the  streamline  region". 

Carman  concluded  that  pressure  drop  through  a  variety  of  packing  materials 

could  be  represented  by 

This  equation  was  formulated  from  information  concerned  with  flow  that  was 
iiighly  1;     ,  transitional,  or  highly  turbulent.  Figure  6  illustrates  how 
this  equation  represented  the  information  available  at  that  time,  the  coordi- 
nates are  identical  to  those  used  by  Burke  and  Plummer. 

Ccr;.ian  personally  conducted  a  test  to  ascertain  the  true  length  of  flow 
path  in  a  bed  of  spheres.  Unaided  visual  observations  showed  that  the  fluid 
path  was  inclined  anywhere  from  0°  to  90°  with  respect  to  the  container  wall, 
and  sometimes  the  fluid  would  follow  an  helical  path.  The  mean  inclination  to 
the  wall  was  concluded  to  be  l»5°« 

A  different  approach  for  estimating  the  hydraulic  diameter  of  a  packed  bed 
was  being  considered  in  the  case  of  packing  materials  such  as  sands  or  crushed 


IOCH 
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ldl9    _^iP£ere^ 


0|  J    Reynold's    ,q    number      [^  ,000 

Fig.  6.     Conclusions  by  Carman  (£)• 

stone,     These  materials  frequently  possess  an  inae terminate  surface  area 
and  the  terra  Ds,   as  defined  in  Table  2,  was  found  satisfactory  for  many 
empirical  correlations.     The  investigation  by  Meyer  and  Vork  (?)   typifies 
this  approach.     They  found  that 


AP     _     /':75Q/<Uo2  (0.67  Vt  -  V) 
*t  gDv?vt 

accomodated  their  crushed  stone  very  well.     Ity  was  expressed  in  terms  of 

Ds  and  V/V-t>   thus  accommodating  wall  surface  and  factors  due  to  variation 

of  free  space.     Use  of  Dg  was  supported  by  analogy  between  diameter  of  a 

packing  unit  and  its  surface  according  to  Sp/Vp  *  6/D^  for  spheres. 

Sullivan  and  Hertel  (lii)   elaborated  on  the  reasoning  that  Carman  used 

to  predict   <S   for  laminar  flow  through  beds  packed  with  spheres  and  textile 

fibers.     They  assumed  that    o  for  packed  beds  should  originate  from  a  basic 

value  of   o  ■  3  for  an  arbitrary  duct  which  existed  within  the  packed  bed. 

They  assumed  that  the  effect  of  "rn  could  be  expressed  by  the  mean  angle  of 

orientation  of  surfaces  with  respect  to  the  overall  direction  of  flow.     Jhey 

expressed  this  mathematically  as    u  ■  3/(sin  6)av.     Their  reasons  for  choosing 
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a  basic  value  of  <&  ■  3  can  be  understood  light  of  the  fsct 

that  the  agreement  with  experimc ntal  results  was  excellent.  An  earlier 
paper  by  Fowler  and  llcrtel  (6)  showed  that  a  basic  value  of  o   "  3  exists 
only  in  the  extreme  case  of  an  infinitely  wide  rectangular  duct.  Sections 
that  might  more  logically  be  present  in  packed  beds  are  the  triangle  or 
square.  Both  of  these  possess  values  for  <$  which  are  much  less  than  3« 
It  is  also  true  that  the  wide  rectangular  duct  suffers  a  large  depression 
in  q  if  the  walls  converge  as  often  as  once  every  ten  times  the  distance 
between  them. 

Table  3.  Values  of  o  f°r  laminar  flow  throuph  earpty  ducts  (7). 


Shape  of  cross-section 


Circular  2.00 

Elliptical 

major  axis  ■  2  x  minor  axis  2.10 

major  axis  *  10  x  minor  axis  2.U2 

Rectangular 

square  1.78 

length  *  2  x  width  l.^ij 

length  ■  10  x  width  2 .65 

infinitely  wide  3«00 

Triangular 

eouilateral  triangle  1.67 

infinitely  high  1.50* 


Included  by  present  author. 


The  experimental  results  of  Sullivan  and  Hertel  are  believed  to  be  the 
most  reliable  in  the  literature.  They  payed  extreme  attention  to  detail 
and  found  that 

££  -  I».50(l  ♦  2  St/3  Sp)2  A  U°Sp  Vt  *  0.5$  percent 
14  *?3 


Ill 

for  their  tests  with  a  few  beds  of  spheres.  Although  this  equation  was 
derived  from  very   general  considerations,  it  is  applicable  only  to  those 
few  beds  of  spheres.  The  wall  correction  factor,  (1  +  2St/3Sp),  was 
arrived  at  by  considering  mathematical  consistency  between  the  circular 
column,  &  m   2,  and  tne  value  for  o  •  u.£0.  This  correction  has  been  found 
slightly  too  large  to  apply  to  greater  values  of  S^/Sp.  Similar  reasoning, 
based  upon  practical  consideration  of  a  great  variety  of  packing  materials, 
has  led  to  another  value.  Sullivan  end  tfertel  maintained  ^UqV^/^  S  <  O.Olli 
thus  insuring  that  totally  insignificant  transitional  effects  were  encountered. 

Wall  effects  as  a  function  of  Hp/D^  has  been  reviewed  by  Perry  (1?).  I  ore 
recent  developments  along  this  line  have  been  accomplished  in  a  series  of 
articles  by  Leva  and  Grummer  (8). 

Oman  and  t'atson  (11)  tested  several  beds  of  different  packing  materials. 
Flow  was  almost  entirely  turbulent.  They  introduced  an  accurate  picture  of 
the  effect  of  "loose  pack"  and  "close  pack"  upon  pressure  drop.  It  was  found 
that  by  plotting  friction  factor  ■  —  — "  ■ vs.  Reynold's  number 

■  on  logarithmic  coordinates  a  best  mean  line  representation  was 

MSp 
obtained.  A   i actor,  (fd/f-^),  was  used  to  allow  for  the  effects  encountered  due 

to  the  differences  in  void  space  offered  by  the  two  methods  of  packing.  One 
important  point  is  that  free  volume  was  included  as  being  raised  to  the  1.7 
power  rather  than  the  3»0  power,  a  direct  contridiction  to  mo*st  previous  con- 
clusions. This  is  in  good  agreement  with  results  obtained  in  the  present 
investigation. 

liorcom  (10)  developed  sn  equation  for  predicting  pressure  drop  in  laminar, 
transition,  or  turbulent  flow  regions.  His  equation  agreed  well  with  experi- 
mental data  for  many  individual  beds  and  was  of  the  form 
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Lt  g  g 

More  detail  is  purposely  excluded  because  reference  was  made  to  some  very 
illusive  terms  such  as  "normal  voids".     The  most  important  consideration  is 
that  this  expression  accommodated  the  various  regions  of  flow  well  for  any 
one  bed.     The  packing  materials  used  were  poorly  defined  so  that  results  of 
his  individual  tests  were  not  useful  in  a  detailed  analysis  of  flow.     It  is 
noteworthy  that  equation  (5)   is  in  general  agreement  with  the  conclusions  of 
Carman  except  for  absence  of  the  term,   (/*§  Jp  U0Vt)       •    An  inspection  of 
Fig.  6  shows  that  pressure  drop  through  spheres  would  have  been  better 
accommodated  had  this  term  been  absent  in  Carman's  equation. 

Illustrations  in  a  text  by  Rouse  (13)  typify  how  the  results  of  studies 
of  flow  about  suspended  objects  tend  to  verify  equation  ($). 

Brownell  and  Katz  (2)  introduced  several  new  concepts  originating  from  a 
comparison  of  pressure  drop  in  packed  beds  to  pressure  loss  in  conduits.     Of 
primary  importance  is  the  discussion  of  the  effect  of  surface  roughness.     Fig- 
ure 7  shows  several  curves  that  were  illustrated;  most  of  these  proposed  curves 
are  approximated  by  equation  (f>)»     The  different  packed  bed  variables  con- 
sidered by  Brownell  and  Katz  are  summarized  as  follows: 

(V-t/V)131  ■  Reynold's  number  function 

(Vt/V)n  -  friction  factor  function 

Sg/sp  *  sphericity 

e/Dp  ■  relative  surface  roughness 

Dp^UoAx  *  modified  Reynold^  number 

(2gDpAP/Lt^U02)  -  modified  friction  factor 


16 

By  means  of  the  Reynold's  number  function  and  the  friction  factor  function, 
a  representation  of  equation  (1)  was  moved  onto  the  pipe  friction  factor 
curve •  Consideration  of  both  sphericity  and  surface  roughness  resulted  in 
satisfactory  representation  for  a  large  series  of  particles  of  "primary 
configuration".  £ome  difficulty  was  encountered  when  "splined"  rings,  or 
particles  of  "second  configuration",  were  considered.  In  the  written  dis- 
cussion to  the  authors,  C.  E#  Lapple  expressed  doubt  as  to  the  possibility 
that  surface  roughness  actually  contributes  to  pressure  drop  in  packed  beds. 
In  answer,  the  authors  agreed  that  surface  roughness  was  of  little  importance. 


friction 
factor 


curvature    suggested   by 
equation  (5) 


s^o// 


10    fold 


Reynold's      number 


Fig.  7.     .Effect  of  surface  roughness  according 
to  Brownell  and  Katz  (2). 
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The  proposals  of  Brownell  and  Katz  were  elaborated  by  Brownell, 
Dombrowski  and  Dickey  (3) •  They  proposed  Figs.  8  and  9  to  represent  the 
Reynold's  number  function  and  the  friction  factor  function.  They  arrived 
at  no  definite  conclusions  concerning  the  effect  of  surface  roughness . 
Intensive  tests  which  were  conducted  on  a  few  beds  tend  to  verify  equation 

(30. 

It  is  of  some  interest  to  note  that  the  Reynold's  number  function  of 
Fig.  8  will  not  lend  itself  to  the  wire  rings  tested  by  the  present  author. 
Figure  9   only  vaguely  suggests  a  friction  factor  function  for  them.  The 
wire  rings  possessed  a  sphericity  of  0.1*2  and  beds  of  them  possessed  82  to 
8JU  percent  void  space. 

MATERIAL  AND  MEOEOBB 

Pressure  drops  through  several  packed  beds  was  observed  in  order  to 
gain  new  knowledge  as  to  the  effect  of  tower  surface  and  packing  density 
upon  pressure  drop.  Pressure  drop  through  less  conventional  packing  materials 
was  also  sought. 

The  Packing  Materials 

Seven  packing  materials,  including  four  different  types  of  packing 
units  were  tested.  The  volume,  surface  area,  and  perimeter  of  each  unit 
was  determined  according  to  a  method  which  was  considered  to  be  direct  and 
accurate . 

The  following  outline  illustrates  the  basic  measurements,  accuracy  of 
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Fig.  8.     Brownell'3  (3)  Reynold » n  number  function. 
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Fig.  9.    Brovmell's  (3)  friction  factor  function. 
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determination,  and  the  derived  information. 
Packing  #1,  wire  ring 

Diameter  of  wire  by  micrometer,  23  arbitrarily 

selected  units:  0.06765^  0.00029  inches 
Volume  by  bouyancy  in  water,  1190  arbitrarily 

selected  units:  3»83  x  10~6  cu.  ft. 
Volume  by  water  displacement,  25700  arbitrarily 

selected  units:  3 #95  x  10""6  cu.  ft. 
Volume  by  water  displacement,  11750  arbitrarily 

selected  units:  h»0h   x  10""6  cu.  ft. 

Final  values: 

volume  -  lj.00  x  10"6  cu.  ft. 
surface  »  2.89  x  10~3  sq.  ft. 
perimeter  ■  0.338  ft.  ■  2L  ♦  nD 

Packing  #2,  glass  ball 

Volume  by  bouyancy  in  water,  109  arbitrarily 
selected  units:  1.185  x  10"^  cu.  ft. 

Volume  by  water  displacement,  £00  arbitrarily 
selected  units:  1.506  x  10~&  cu.  ft. 

Final  values: 

volume    -  1.502  x  10"**  cu.  ft. 
surface   •  6  v/D,  D  -  0.0660  ft. 
perimeter  «  0.207li  ft.  -  ffD 

Packing  #3>  clay  Berl  saddle 

Width  by  machinist's  rule,  20  arbitrarily  selected 

units:  0.085ii  ft. 
Width  by  steel  tape,  30  arbitrarily  selected  units 

placed  edge  to  edge:  0.08ii7  ft. 
Volume  by  bouyancy  in  water,  96   arbitrarily 

selected  units:  l.Uili  x  10""^  cu.  ft. 
Volume  by  water  displacement,  275  arbitrarily 

selected  units:  1.37ii  x  10~^  cu.  ft. 
Volume  by  water  displacement,  500  arbitrarily 

selected  units:  l.ii68  x  10""*1  cu.  ft. 

Final  values: 

volume    •  l.iiU5  x  lcH*  cu.  ft. 

surface   *  0.03^3  sq.  ft.  (from  mfg.  data*) 

perimeter  *  0.533  ft.  -  2itD 


Manufactured  by  the  Maurice  A.  Knight  Company,  Akron,  Ohio, 
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Packing  #h,   clay  Berl  saddle 

!th  by  steel  tape,  itO  arbitrarily  selected  units 
placed  edge  to  edge:     O.OIj27  ft. 
Volume  by  water  displacement,   21)00  arbitrarily 
selected  units:     2.5k  x  10-S  Cu.  ft. 

Final  values j 

volume         ■  2.51  x  lcH*  cu.  ft. 

surface       *  8.80  x  10"^  sc.  ft.  (from  mfg.  data*) 

perimeter  -  0.268  ft.  •  2itD 

Packing  #5*  clay  Raschig  ring 

Volume  by  water  displacement,  330  arbitrarily 

selected  units?     2.68  x  1CT^  cu.  ft. 
Diameter  by  steel  tape,   58  arbitrarily  selected 

units  placed  side  by  side*     0.0860  ft, 
Length  by  steel  tape,  95  arbitrarily  selected 

units  placed  end  to  end:     0.0873  ft* 

Final  value st 

volume         *  2.68  x  10"**  cu.  ft. 

surface       ■  O.Oii59  sq.  ft. 

perimeter  -  0.630  ft.  ■  2  irD  ♦  2L  -  2  wt 

Packing  £6,  clay  Raschig  ring 

Volume  by  water  displacement,   2000  arbitrarily 

selected  units t     lull  x  10"?  cu.  ft. 
Diameter  by  steel  tape,  °1»  arbitrarily  selected 

units  placed  side  by  side:     O.Oli35  ft. 
Length  by  steel  tape,  16?  arbitrarily  selected 

units  placed  end  to  end:     O.Olilili  ft. 

Final  values: 

volume         ■  ij.ll  x  10""?  cu.  ft. 

surface      -  0.01165  sq.  ft. 

perimeter  -  0.312  ft.  -  2  irD  +  2L  -  2  trt 


^Manufactured  by  the  Maurice  A.  Knight  Company,   Akron,   Ohio. 
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Packing  f"7,  metal  Raschig  ring 

Diameter  by  micrometer,   25  arbitrarily  selected 

units:     1.0115*  O.OOI49  inch 
Length  by  micrometer,   25  arbitrarily  selected 

units:     1.0031**  0.0005  inch 
Density  by  Testphal  balance,  water  displacement, 

5  arbitrarily  selected  units:     7.876I4-  0.0135  go. 

per  cu.  cm. 
Mass  of  total  supply,  698  -  81*56—  6  grn. 

Final  values: 

percent  void  "  92* 

surface     ■  62.7  sq.  ft.  per  cu.  ft.* 

perimeter    -  0.681  ft.  ■  2  i»D  ♦  2L  -  2  irt 

The  calculations  for  perimeter  are  indicated  primarily  to  illustrate 
how  the  perimeter  is  defined  for  packing  materials.  Determination  of  the 
perimeter  of  the  Berl  saddle  is  considered  only  an  approximate  method;  all 
others  are  exact  according  to  the  present  definition.  Auxiliary  measure- 
ments, such  as  thickness  of  the  Raschig  rings  or  length  of  the  wire  rings, 
were  observed  to  agree  with  the  above  conclusions. 

The  Flow  System 

Three  steel  pipes  of  differing  diameter  were  used  as  columns.  E?ch 
of  these  was  thirty-six  inches  in  length.  The  packed  zone  included  the 
entire  length  of  each  pipe  while  pressure  drop  measurements  were  taken 
from  pressure  taps  located  twenty-four  inches  apart  and  six  inches  from  the 
ends. 

Plate  II  illustrates  the  exact  flow  system  used  and  the  location  of 
the  various  metering  instruments. 


*Manufacturer»s  information:  Metallo  Gasket  Company. 
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The  pressure  taps  in  the  four  inch  column  were  different  from  the 
piezometer  rings  usee!  in  the  three  and  six  inch  columns.  These  were  so 
constructed  to  facilitate  a  later  study  of  counter-current  flow.  No 
difference  in  results  was  noted  that  could  be  attributed  to  the  difference 
in  style  of  the  pressure  taps. 

The  diameters  of  the  columns  were  determined  by  filling  each  column 

with  water  and  noting  the  amount  required  for  the  space  between  the  pressure 

taps.  Column  diameters  and  sectional  areas  thus  determined  were  as  follows: 

Column        Diameter^  inches        Area,  square  feet 

three  inch         3.10  0,05207 

four  inch  h.06  O.O898 

six  inch  6. 08  0.2019 

The  porosity  of  each  packed  column  was  determined  by  the  sane  process 
used  to  determine  column  diameter.  The  porosity  of  the  bed  of  metal  rings 
was  the  exception,  its  void  fraction  was  ascertained  from  information  pub- 
lished by  the  manufacturer. 

Fluids  and  Flow  Measurement 

Three  fluids  were  used:  SJV.E.  ^60  oil,  water,  and  air. 

The  density  of  the  S.A.E.  /'60  oil  was  determined  by  T.'estphal  balance. 
The  balance  was  calibrated  against  water  samples  at  various  temperatures. 
The  oil  samples  of  varying  temperature  were  then  tested.  The  balance  was 
found  to  be  very  sluggish  when  measuring  the  density  of  the  oil;  this  was 
overcome  by  allowing  sufficient  time  for  the  balance  to  react.  No  attempt 
was  made  to  control  temperature  closely,  temperature  being  read  on  the  plumet 
of  the  balance,  since  it  was  not  anticipated  that  serious  effects  of  convec- 
tion would  be  existant  at  the  range  of  temperatures  encountered.  Results 
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were  reproducable  to  -within  ±0.0005  gram  per  cubic  centimeter.  The  results 
used  are  as  follows* 


..   o.  °C 

Density 3   lb,  j     •  cu.  ft. 

20 

56.28 

25 

56.11 

30 

55.95 

3* 

55.75 

Lo 

55.56 

Viscosity  of  the  oil  was  determined  by  the  Kansas  State  Highway  Department 

at  different  temperatures.  Repeated  checks  on  samples  that  contained 

possible  impurities,  such  as  sludge,  emulsified  water,  or  emulsified  air, 

showed  that  little  error  resu'ted  from  the  presence  of  these  impurities. 

Results  of  these  tests  are  as  follows? 

Date     Temp.,   Viscosity,    Possible 
_      °F     centistokes   impurities 

12-16-50   70  1636.6  dissolved  water 

80  1038.5       "      » 

100  h$h*72                n              " 

210  ?5.5l      "     " 

l-ij-51    100  W»5.03  sludge  and  water 

100  U33.68      «   »'   ■ 

2-6-51    100  131.6       ■   "  air 

100  1)32 .2 


it   tt   n 


The  impurities  are  noted  to  have  affected  viscosity  very  little  over  a  period 

of  two  months  usage.  The  average  viscosity  of  samples  at  100°  F  for  tests 

on  12-16-50  and  l-h-51  were  used  as  a  basis  for  calculations.  The  trend  of 

viscosity  with  temperature  was  determined  from  results  of  12-16-50;  the 

logarithm  of  absolute  viscosity  was  found  to  vary  linearly  with  the  inverse 

cube  of  absolute  temperature.  The  validity  of  the  relationship  is  illustrated 

thust 

Temp,  range,   °F  70-80       80-100       100-210 

Alog^/  A(looo/°R)3       0.5^3     0.5195     0.531- 
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The  following  viscosity  information,  being  between  80  and  110°  F,  was 

derived  from  the  slope  determined  for  80  to  100°  F  of  0.261  lb.  per  ft  .-sec. 

Temp.  °C  30  32  3)4  36  38  10  b2 

Vise.  rel. 

to  100°  F      1,771*    3  .$20    1.310    1.132    0.985    0.859   0.750 

Viae,   lb. 

per  ft-sec.  0.163        0.397        0.31*2        0.295        0.257        0.221*      0.196 

These  viscosities  were  used  for  all  flow  calculations.     Flow  of  the  oil  was 

measured  by  time  required  for  a  weighed  quantity  of  oil  to  flov;  from  the  system. 

The  density  and  viscosity  of  the  water  were  taken  from  the  Handbook  of 
Chemistry  and  Physics  (8).     The  -water  used  was  obtained  directly  from  the 
Manhattan  city  supply. 

The  density  and  viscosity  of  dry  air  were  obtained  from  the  same   source 
as  the  information  for  water.     Corrections  for  moist  air  were  applied  as 
follows: 

dens,  moist  air  *  dry  air  (1  -  0,6l  abs.  hum.) 

vise,  moist  air  *  vise,  dry  air'  -  abs.  hum.  (vise,  dry  air  -  vise. 

wat.  vap.) 

The  estimation  for  density  of  moist  air  is  an  approximation  which  is  good 

for  low  values  of  absolute  humidity.  No  mumidities  over  one  percent  were 

encountered.  The  estimation  for  the  viscosity  of  moist  air  is  based  on  the 

assumption  that  viscosities  are  additive  with  respect  to  weight  percent.  The 

calculated  viscosity  was  never  less  than  99  1/2  percent  of  the  viscosity  of 

dry  air.  The  air  was  obtained  from  the  compressed  air  supply  of  the  Kansas 

State  College  of  Agriculture  and  Applied  Science.  Humidity  was  measured 

with  a  sling  psychrometer  and  interpreted  according  to  the  psychrometric 

chart  from  the  textbook  of  Badger  and  .\5cCabe  (17)  • 
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How  of  air  and  water  was  measured  with  a  flow  nozzle  made  by  expanding 
one  end  of  a  short  length  of  brass  pipe.     The  nozzle  was  0,0552  ft,  in  dia- 
meter and  was  mounted  in  a  one  inch  steel  pipe  line.     Impact  and  static 

pressure  taps  were  located  at  the  exit  of  the  nozzle.     The  nozzle  coeffi- 

1/2 

cient  was  found  to  be  constant  at  W/ /°  (&li)      '  ■  0.00513  for  the  range  of 

flow  studieu.     Calibration  was  maae  with  water  where  7»  ■  flow,  lb,  per  sec; 
P    =*  density,  lb,  ^or  cu.  ft,;   AH"  head  loss,   inches  of  fluid.     This 
meter  was  later  calibrated  with  air.     The  coefficients  were  found  to  agree 
within  four  percent  for  the  different  fluids.     The  above  mentioned  coeffi- 
cient corresponds  to  a  discharge  coefficient  of  £.bout  93  percent.     The  flow 
nozzle  used  is  seen  in  place  in  Plate  I,   together  with  other  nozzles  of 
similar  construction. 

Manometry 

Manometers  were  used  to  measure  pressure  differential  caused  by  flew 
through  the  nozzles,  pressure  drop,  and  pressure.  Inverted  manometers 
were  used  when  water  or  oil  was  in  the  system;  water  filled  manometers  were 
used  for  air  flow  and  pressure  drop,  ilcrcury  filled  aanometers  were  used  to 
measure  pressure.  All  manometers  were  the  "U"  tube  type  and  were  calibrated 
in  inches  of  fluid  displacement.  They  are  illustrated  in  Plate  II, 

Each  manometric  reading  was  interpreted  so  as  to  include  the  secondary 
effects  of  air  as  a  second  fluid  and  the  difference  between  the  density  of 
the  fluid  in  the  system  and  the  density  of  the  fluid  in  the  manometer.  The 
latter  effect  was  significant  when  the  temperature  of  the  oil  approached  1»0 
degrees  centigrade. 
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Thermometry 

Thermometers  were  located  as  indicated  in  Plate  II.  The  column  tempera- 
ture was  determined  as  the  median  temperature  between  inlet  and  outlet 
points.  The  temperature  change  did  not  exceed  two  degrees  centigrade  for 
any  one  run.  The  temperature  of  the  flow  meter  was  assumed  to  be  the  same 
as  the  temperature  indicated  by  the  inlet  thermometer.  All  thermometers 
were  checked  against  a  precision  thermometer  so  that  any  one  temperature 
reading  could  be  considered  accurate  within  ^0.1  to  -0,2  degrees  centigrade* 
The  possibility  of  a  temperature  gradient  at  right  angles  to  the  flow  path 
was  considered.  It  was  found  that  a  gradient  of  only  one  degree  centigrade 
existed  when  oil  in  the  reservoir  was  at  a  temperature  of  35  degrees  centi- 
grade. Since  oil  moved  far  more  rapidly  through  the  system  than  it  did 
through  the  reservoir,  it  is  assumed  that  no  measurable  gradient  existed  any- 
where within  the  flow  system. 

Sample  Experimental  Data 

When  oil  was  used  in  the  system,  circulation  was  maintained  for  one 
half  hour  at  each  different  rate  to  insure  equilibrium  in  the  manometers 
and  to  insure  thermal  equilibrium  in  the  system.  Two  hours  were  allowed 
for  initial  equilibrium  for  each  series  of  runs.  Apparent  equilibrium  was 
reached  in  half  of  the  allowed  times.  The  following  information  was 
gathered  twice  in  succession,  sometimes  three  times,  in  order  to  determine 
pressure  drop,  flow  rate,  viscosity,  and  density  for  each  single  "run". 
Since  the  time  pattern  for  each  reading  was  symmetrical  about  the  flow 
reading,  a  direct  numerical  average  of  results  was  made. 
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Run  #87*     (first  half) 

manometer  temp.  •  23»7°  C 

lower  manometer  leg  »  lli.5>5>  inches 
upper  manometer  leg  •  39 #72  inches 


outlet  temp.  ■  37.2°  C. 

flow  weight  -  1920  less  371  gm 

flow  time  ■  2iul  sec. 

time  -  08x36  3/2 

inlet  temp.  ■  38.2°  C 

upper  manometer  leg  ■  39»69  inches 
lower  manometer  leg  "  llj.52  inches 


f  20  sec. 

}  nearly  simultaneous 
]  20  sec. 


simultaneously 

\     20  sec. 

)     nearly  simultaneous 
1     20  sec. 


manometer  temp.  »  2l4.0°  C 

Most  of  the  runs  were  made  during  winter  months  in  a  large  laboratory 
that  opened  out  of  doors,  and  it  was  not  uncommon  for  the  opening  of  a  door 
to  cause  room  temperature  to  suddenly  drop  one  half  to  one  degree  centigrade. 
This  would  cause  the  air  in  the  inverted  manometer  to  contract  and  draw  both 
legs  up  a  fraction  of  an  inch.     The  differential  readings  did  not  vary  by 
more  than  0.0f>  inch  in  any  case.     Runs  identified  by  alphabetical  symbols 
were  accomplished  less  systematically;  they  were  noted  to  yield  the  same 
results  as  the  remainder  of  the  tests. 

Runs  with  water  were  similar  to  those  with  oil  except  that  only  five 
or  ten  minutes  were  required  for  equilibrium  and  that  flow  was  measured  with 
a  nozzle.     An  entire  series  of  readings  could  be  taken  within  30  seconds, 
thus  a  time  schedule  was  not  maintained.     Double  readings,   as  below,  were 
usually  taken.     Sometimes  a  fine  oil  ring  in  the  manometer  tube  facilitated 
reading  so  well  that  one  reading  was  considered  sufficient. 

Run  'r79: 

col.  man.  temp.  »  21.3°  C 

lower  leg  column  man,  ■  I8.I45  inches 

upper  leg  column  man.  »  29.20  inches 
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Run  #79  (cont.) 

outlet  temp.  ■  26.2°  C 

inlet  temp.  -  26.3°  C 

flow  man.  temp.  »  21.3°  C 

static  flow  leg  3  -11.3  inches 

impact  flow  leg  •  -*1$#35>  inches,   time  ■  13:2b 

impact  flow  leg  "  +15.35  inches 

static  flow  leg  ■  -11.35  inches 

inlet  temp.  -  26.1°  C 

outlet  temp.  -  26.3°  C 

upper  leg  column  man.  •  29*00  inches 

lower  leg  column  man.  *  l8.1jj>  inches,   time  *  13 :2£ 

Runs  with  air  required  more  information  that  runs  with  water*     Equili- 
brium was  reached  so  rapidly  that  a  time  schedule  was  not  considered  useful. 
Duplicate  readings  were  made  for  each  run  as  illustrated  below.     Some  fluctua- 
tions in  readings  were  noted,  but  they  were  so  rapid  that  a  time  schedule 
for  making  readings  would  not  have  been  capable  of  capturing  the  average 
reading  any  better  than  a  fast  scanning  of  all  instruments. 
Run  #105:     (first  half) 

impact  leg  flow  man.  ■  -»-l.i<2  inches 

static  leg  flow  man.  »  ♦2.83  inches 

flow  man.  temp.  -  22.8°  C 

meter  side  of  gage  pressure  man.  for  static  meter  tap  »  -0.30  inches 

atmospheric  side  of  gage  pressure  man.  for  static  meter  tap  «  +0.16  ins. 

inlet  temp.  -  23 .2°   C 

outlet  temp.  ■  23.0°  C 

upper  leg  col.  man.  ■  -0.$9  inches 

lower  leg  col.  man.  *  +0.81  inches 

col.  man.  temp.  ■  23 »Q°  C 

system  side  of  gage  pressure  man.  for  upper  column  tap  ■  +0,83  inches 

atmospheric  side  of  gage  pressure  man,  for  upper  column  tap  »  +1.09  ins. 

time  •  16 j 12 

Barometric  pressure  and  humidity  of  the  exit  air  supplemented  this  information. 

Fluctuations  in  manometer  readings  never  produced  discrepancies  greater  than 

2  percent  of  the  manometer  displacement  for  any  one  run. 
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Schedule  for  Each  Bed 

The  beds  were  packed  by  introducing  about  five  to  ten  percent  of  the 
required  packing  material,  settling  this  by  rapping  the  column,  and  then 
introducing  another  five  to  ten  percent  of  the  required  packing  material. 
The  mixture  in  bed  .fllj  was  introduced  in  individual  portions  that  repre- 
sented the  simplest  subdivision  of  the  mixture.  Free  space  was  measured 
immediately  after  packing  each  column,  and  after  all  runs  were  completed* 
No  settling  of  the  packing  during  runs  was  noted,  "Water  was  the  first 
fluid  used,  air  was  next,  then  oil. 

The  entire  system  was  flushed  with  carbon  tetrachloride  and  dried  after 
tests  with  oil  were  completed.  Then  the  system  was  flushed  with  water. 
Some  oil  remained  in  the  system,  but  never  any  more  than  enough  to  produce 
a  thin  oil  slick  on  top  of  the  water, 

THEORETICAL  DEVELOPMENTS 

General  Considerations 

The  first  consideration  was  that  of  locating  a  more  useful  equation 
than  that  resulting  from  dimensional  analysis.  It  was  decided  that  Korcomfs 
(10)  representation  should  adequately  determine  the  relationship  between 
pressure  drop  and  fluid  variables. 

A-£-  k^£  +  K  {li!2L  (5) 

H      g  g 

Equation  (f>)  has  properties  such  that  laminar  and  turbulent  flow  may 
be  scrutinized  independently.  Actually,  few  experimental  data  concerning 


EXPLANATION  OF  PLATE  II 

Sketch  of  the  flow  system  showing  instrument  location 

C-l  The  six  inch  column. 

C-2  The  three  inch  column. 

C-li  The  four  inch  column. 

E  Exit  manifold  for  all  three  columns. 

F  Location  of  the  flow  meter  used  to  measure  the  flow  rates  of  water  and  air. 

I  Inlet  manifold  for  all  three  columns. 

M-l  Inverted  manometer  used  for  measuring  the  rate  of  water  flow. 

14-2  Mercury  filled  manometer  used  to  measure  pressure  in  the  columns. 

14-3  Water  filled  manometer  used  to  measure  pressure  drop  for  runs  with  air. 

i  -h  Inverted  manometer  used  to  measure  pressure  drop  for  runs  with  oil  or  water. 

M-5  "Water  filled  manometer  used  for  measuring  the  rate  of  flow  of  air. 

ii-6  Mercury  filled  manometer  used  for  measuring  air  pressure  in  the  flow  meter. 

P  Piezometer  rings  and  pressure  taps. 

R  Reservoir  for  fluid  being  circulated. 

S    Positive  displacement  pump,  eccentric  gear  type,  used  to  circulate  water 
or  oil. 

T    Thermometers  used  to  measure  the  temperature  of  the  inlet  and  outlet 
streams. 

-/    Portions  of  the  pipe  system  that  were  closed  to  the  circulating  fluid. 


PLATE  II 


EXPLANATION  OF  PLATE  III 
Photograph  of  the  flow  system. 
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PLATE  III 


EXPLANATION  OF  PLATE  IV 
Photographs  of  the  packing  units. 

Fig.  10.  Upper  left,  one  half  inch  clay  Berl  saddle. 

Upper  right,  one  inch  clay  Raschig  ring. 

Center  left,  -wire  ring. 

Center  right,  glass  ball. 

Lower  left,  one  half  inch  clay  Raschig  ring. 

Lower  right,  one  inch  clay  Berl  saddle. 
Fig.  11.     The  mixture  tested  in  bed  #llu 
Fig.  12.     The  raetal  Raschig  rings* 


PLATE     IV 


35 


<      <       4* 

^^^^    ^^^^     ^^^^^ 

*       *'     # 

**^? 

o     o     o 

*  *  * 

o    o    o 

&  &  <* 

o    o    o 

<9   <9    <9 

*  *  * 

*  «r* 

*.  *  * 

Jt  *  ** 

/ 

***: 

i 

\*  *  * 

Fig.     1 1 


Fig.       10 


Fig.    12 


36 

truly  turbulent  flow  -were  found.  This  required  that  laminar  flew  be  studied 
first.  Development  of  a  reliable  prediction  for  an  equivalent  of  MkM  made 
possible  direct  evaluation  of  an  equivalent  for  "Kw  from  data  obtained  for  some- 
what transitional  flew. 

Laminar  Flow*  In  order  to  obtain  agreement  with  equation  (1),  which 
summarizes  dimensional  analysis,  it  was  necessary  to  approximate  certain  bed 
variables. 

Since  fluid  and  empty  column  variables,  A  p,  h^tylM »  P  >     Uo>  were  known,  a 
limited  number  of  others  was  required.  Carman's  (5>)  representation  for  wmn 
was  adopted  to  replace  "D"  and  was  later  found  to  be  the  proper  substitution. 
This  was  defined  by  equation 

m  -  V/Sp  (2) 

On  first  consideration,  the  classic  approximation  for  area  of  flow,  used  by 
Blake  (l)  and  many  other  investigators  was  thought  to  be  useful.  This 
amounted  to  reducing  the  column  cross- section  by  the  void  fraction.  Later 
considerations  found  this  quite  valueless. 

The  final  list  of  bed  variables,  in  addition  to  the  approximation  for 

hydraulic  radius  was  concluded  to  be: 

a  »  effective  area  of  flow 
column  cross-section 

r  »  e*f  ecfriVG  length  of  flow  path 
length  of  the  bed 

«  effective  hydraulic  radius 
estimated  hydraulic  radius 
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0\  ■  constant  derived  from  geometrical  nature  of  the  flow  path 
Sfc/Sp  ■  ratio  between  column  surface  and  packing  surface 
Pressure  drop  in  laminar  flour  has  long  been  known  to  be  independent  of  sur- 
face roughness. 

Turbulent  Flow,  The  same  terms  that  were  related  to  laminar  flow, 
except  for  o,  »  were  considered  to  be  applicable  to  turbulent  flow.  This 
may  not  have  been  an  exact  assumption  because  for  instance,  the  effective 
area  of  flow  available  to  a  turbulent  stream  might  be  different  from  the 
area  available  to  a  stream  in  laminar  motion. 

The  effect  of  surface  roughness  was  expected  to  become  evident  in  the 
turbulent  region.  Additional  terms  to  be  considered  included: 

6xm   constant  derived  from  geometrical  nature  of  the  flow  path 

e/m  ■  height  of  surface  protrusion 
estimated  hydraulic  radius 

The  Basic  Equation,  Restatement  of  equation  (Jj>)  produced  equation 

(6). 

2 


A? 


&*£ .  *  */ftK£J&! .  ±i/h  s)       (6) 

'  gm^  az?  'l^j      >  aZss     I  Sp    m) 


tf   was  used  to  represent  the  arbitrary  functions.  Examination  of  equation 
(6)  showed  that  attempts  to  solve  for  all  of  the  bed  variables  from  pressure 
drop  information  would  be  futile,  A  simplier  form  was  adopted  for  further 
analysis  of  pressure  drop, 

4£  -  A»  ^  Uo  +  n*/°uo2  (7) 

Both  A1  and  B*  could  then  be  quickly  evaluated  for  each  bed.  It  was 
felt  that  any  randomly  packed  bed  should  possess  values  for  o,  and  03.  that 
depended  on  the  same  bed  properties  that  would  determine  a,  r,  and  2,  therefore, 
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all  five  of  these  terms  would  be  accountable  to  some  single  bed  variable. 

The  v.idth  of   Packing  Unit,  Satisfactory  determination  of  A«  and  B* 
was  obtained  by  introducing  the  concept  of  packing  width.  The  width  of 
the  packing  is  a  derived  property  which  is  not  necessarily  related  to  the 
nominal  diameter  of  the  packing  unit. 

The  packing  width  was  defined  as  the  surface  area  of  the  packing 
divided  by  the  perimeter  representing  boundries  which  must  be  circumvented 
by  the  fluid.  This  Yras  expressed  mathematically  as: 

In  order  to  facilitate  a  precise  estimation  for  Cp,  this  definition  was  for- 
mulated: the  packing  perimeter  consists  of  the  locus  of  tangent  points  to 
the  packing  that  would  be  generated  by  a  line  which  moved  throughout  the 
packed  bed  remaining  oriented  parallel  to  the  column  wall. 

Comparison  of  "w"  to  "m"  yielded  a  variable  that  uniquely  measured 
the  degree  to  which  the  fluid  path  would  be  distorted  and  blocked.  Scaled 
illustrations  such  as  Figures  1,  2,  and  3  showed  that  flat  packing  units 
might  produce  a  bed  containing  dead  spaces  to  such  a  degree  that  a  bed  of 
lesser  porosity  comprised  of  somewhat  spherical  units  should  produce  no 
more  pressure  drop.  The  magnitude  of  w/m  seemed  to  parallel  this  effect. 
As  a  result  of  these  observations,  it  was  considered  feasible  that  w/m 
would  index  basic  changes  in  the  bed  structure  so  that  the  three  terms, 
w/m,  St/Sp,  and  e/m  could  completely  describe  a  packed  bed. 

Application  of  the  hypotheses 

Preliminary  Analysis.  The  data  of  several  investigators  was  used  to 
supplement  the  experimental  results  of  this  investigation.  Preliminary 
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considerations  showed  that  w/m  very  decidedly  indexed  A'    and  B*  of  equation 
(7).     A  graphical  representation  similar  to  Plate  III  showed  that  the 
limiting  value  for  A1   as  w/m-*0  was  50/9  •     This  initial  representation  also 
showed  that  B*  varied  directly  as  did  w/ra  for  large  columns. 

The  method  of  correcting  for  wall  effects  that  was  used  by  Sullivan 
and  Hertel  (ll;)   was  adopted.     Thus,  it  was  assumed  that 


*P  -  A  ^LSSj  (1  ♦  CoSt/Sp)     *  B*  l-j£.  (8) 


would  accomodate  wall  effects  for  laminar  flow.     That  this  is  true  is 
illustrated  in  Fig.  1$ .     This  equation  produced  precise  correlation  for  wall 
effects  in  the  case  of  spheres,   and  good  correlation  for  all  packing  materials 
for  ratios  of  3^/Sp  from  0.0198k  to  0.3o£.    B*  was  found  more  nearly  uniformly 
dependent  upon  w/D^  than  upon  S-j/Cp.     The  therm,  B    -#-w/m,  was  found  to  depend 
upon  w/D+  but  no  longer  upon  w/m.     The  final  flow  equation  took  this  form* 

M  -  k    ^.Vo  (1  *  0.6St/Sp)2  *  B  /**•**  (9) 

with  A  depending  on  w/ia  and 
B  depending  on  w/D^ 

Contributions  of  e/m  to  pressure  drop  could  not  be  isolated  by  comparison 

to  values  for  "e"  which  were  published  by  Browneli  and  Xatz  (2).  Thus,  it 

was  assumed  that  normal  roughness  should  not  affect  pressure  drop  through 

randomly  packed  beds. 

Istinations  for  Perimeter.  The  contributions  to  perimeter  offered  by 

many  packing  units,  the  sphere,  the  wire  or  cylinder,  the  Raschig  ring,  the 

prism,  and  the  cube,  were  noted  to  be  independent  of  orientation  within  the 

bed.  The  Berl  saddle  was  noted  to  yield  different  perimeter  with  each 

orientation. 
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Observing  the  Berl  saddle  from  various  directions  showed  that  the 
outer  edges  constituted  the  perimeter  from  some  views  while  part  of  these 
outer  edges  ceased  to  contribute  and  other  elements  of  perimeter  appeared 
in  other  views.  For  this  reason,  the  outer  edges  were  felt  to  approximate 
the  wean  perimeter  of  the  Berl  saddle.  The  perimeter  of  the  'saddle"  tested 
by  Brownell  and  co-workers  (3)  was  approximated  by  assuming  that  the  units 
were  manufactured  from  square  blanks  and  that  the  edges  of  the  square  con- 
stituted an  equivalent  of  the  final  perimeter* 

Table  h   shows  the  exact  method  used  to  estimate  the  perimeter  of  each 
different  packing  unit* 

Table  lw  Perimeter  of  some  packing  units* 

Unit  ;  Perimeter 

Sphere 

Viire  or  cylinder 
Cube 

Hexagonal  prion 
Raschig  ring 
Berl  saddle 
"Saddle"  of  (3) 

Laminar  Mow.  Values  for  A  were  calculated  from  information  for  each 
individual  experimental  run*  For  each  bed  tested*  the  logarithmic  mean 
value  of  A  was  Determined.  Table  5  includes  these  results  together  with 
other  important  information,  be tailed  lists  of  the  calculated  results  are 
included  in  the  appendix* 
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•le  5.   Laminar  flow,  moan  values  ox  A  for  each  bed. 


Ref . 

Packing 

-id 

t/sB 

w/m 

.1. 

A 
.  .1. 

',]> 

Glass   ball 

kl.  2 
93.1 

O.llk 

c.56 

10.21 

10.81 

0.945 
1.310 

Smooth   3addle 

0.041 

5.35 

11.11 

6.46 

tt 

Rou^h  saddle 

93.5 

0.04 

5.02 

7.65 

6.26 

0.950 

ft 

rl   saddle 

72.5 

0.113 

6.33 
6.% 

6.58 

9.16 

0.718 

tt 

Raschig  ring 

70.7 

0.121 

7.69 

9.16 

O.t 

(5) 

Hexagonal   prism 

37.7 

0.066 

9.39 

12.32 

11.66 

1.057 

n 

n "                  it 

4.2.6 

0.065 

7.66 

9.30 

10.17 

0.9l4 

ft 

Cube 

31*4 

0.075 

11.11 

4.10 

13.36 

1.055 

it 

tt 

39.7 

0.076 

9.10 

10.80 

11. 40 

0.947 

0.065 

ft 

it 

1*4.8 

0.074 

740 

.62 

9.97 

Ok) 

Glass   ball 

39.08 

0.01965 

9.353 

11.56 

11.63 

11.64 

0.996 

ft 

it           it 

39.04 

0.01984 

9.366 

11.63 

0.999 

ft 

ft           tt 

39.11 

0.09314-7 
O.09367 

9.341 

ll.6l 

11.62 

0.999 

ft 

tt           ft 

39.2Ji 

9.290 

11.64 

11.57 

l.OOo 

Auth. 

-Vire  ring 

81.6 

0.059 

1.39 

6.14 

6.20 

0.990 

ft 

it          it 

83.6 

0.131 

1.21 

6.27 
5.96 

6.11 

1.026 

n 

it          tt 

63.2 

O.098 

1.25 

6.13 

0.972 

tt 

Glass   ball 

42.3 

0.225 

e.20 

10.72 

10.62 

1.009 

ft 

ft                         Tf 

44.0 

0.305 

7.6k 

9. 60 

10.16 

•  ?65 

tt 

irl   saddle 

77.0 

0.265 

4.54 

9.12 

7.95 

1.147 

tt 

tt              n  . 

71.8 

0.118 

5'.?4 

8.52 

8.86 

0.959 

n 

Raschig  rin^ 

68.9 

0.157 

5.63 

11.63 

1  .67 

1.364 

n 

tt               tt 

92 

0.126 

443 

10.08 

7.86 

1.297 

t» 

fi               ti 

92 

0.126 

443 

9.59 

7.66 

1.217 

n 

ti               »f 

59.5 

0.103 

7.20 

9.74 

9.61 

0.993 

H 

*cure;   wire   ring, 
glass   ball,    Berl 

saddle,    Raschig  r. 

53.0 

0.157 

3.57 

8.39 

7.37 

I.140 
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This  rearranged  form  of  equation  (8)   was  used  for  estimations: 

A  .  *pg*2 B,   „ 

L^uUo   (1  ♦  O.oS^Sp)2 

where  N  -  ^U0  rn^a      and  B'   -  B*/(l  +  0,6S|/Sp)2 
The  values  for  B    which  were  used  in  this  estimation  were  obtained  by  initial 
observation  of  turbulent  flow  data  for  the  bed  being  considered*     Since  B 
was  only  approximate,  the  second  term  on  the  right  was  not  allowed  to  exceed 
5  percent  of  the  estimated  value  for  A* 

Plate  V  shows  how  A  varies  with  w/m«     No  theoretical  considerations  to 
explain  the  linear  relationship  on  semilogarithmic  coordinates  were  deduced. 

The  data  of  Sullivan  and  Hertel  (H4)  and  the  intercept  for  A  ■  5&/9  at 
w/m  »  0  were  used  to  determine  this  empirical  relationship: 

A  .  52  (10)0.03^30w/m 
9 

or  log  A  -  O.7I4I47  ♦  0.03ii30w/m  (10) 

Equation  (10)  was  termed  the  "mean  line"  value  for  A, 

Turbulent  Plow «     Values  for  B     were  solved  by  using  this  modified  form 

of  equation  (8)  : 

B*  -    iZg     ■  -  A'/N  -  f  -  A'/N 
/°0  * 


V     no- 


where N  -  m  UQ  P  /jjl   and  A»  -  A(l  +  O^Sf/S  ) 
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B     was  then  converted  to  B  according  to  B  •  B*  +w/m#     Values  for  A1  were 
obtained  by  use  of  equation  (10)  #     This  was  done,   even  when  the  true  value 
was  known,   for  the  sake  of  maintaining  a  consistent  approach  for  all  infor- 
mation that  was  at  hand. 
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The  term,  A'/N,  wts  allowed  to  exceed  the  calculated  value  for  B     in  a 
few  cases,  but  was  gnerally  restrained  to  a  maximum  value  equal  to  the 
estimated  value  for  B  ,     The  logarithmic  mean  value  for  B  for  each  bed  was 
determined  and  placed  in  Table  6  along  with  other  important  information, 

Plate  VI  shows  how  B  varies  with  w/Dt,,     Definition  of  the  mean  line 
value  for  B  according  to 

B  -  0.2^(10) -1*766w/Dt  (11) 

was  found  quite  satisfactory.     This  was  established  by  first  noting  that  B 
approached  0,2£  as  w/D^  approached  zero,   and  by  determining  the  line  that 
would  pass  through  this  point  presenting  the  closest  approximation  to  the 
mass  of  information, 

r/m  As  a  Unique  Variable,     If  w/m  were  not  a  unique  variable,  A  and  B 
would  show  dependence  on  the  other  variables,     Plate  VII  shews  that  the 
deviations  of  the  constants  from  their  mean  line  values  do  not  arise  from 
either  the  percent  voids  or  S-j/Sp,     The  fact  that  independence  of  S^/S- 
existed  showed  that  Dp/D^  would  not  explain  deviations  eitr.er,  Dp/D-t  is 
approximately  proportional  to  S^/Sp  for  any  given  packing  material. 

The  variable,  w/m  was  thus  established  as  the  sole  criterion  for  A  and 
*/D+  as  the  criterion  for  B,     Deviations  from  the  mean  line  values  were 
attributed  to  normal  experimental  errors  such  as  might  arise  from  insuffi- 
cient column  length  when  large  packing  units  were   tested  in  small  diameter 
columns. 

Friction  Pactor  and  Reynold's  Number,     It  was  desired  to  represent  all 
experimental  information  on  a  single  friction  factor  vs,  Reynold's  number 
plot.     In  order  to  do  this,    certain  considerations  of  the  variability  of  A 


Table  6.  "  Turbulent  flow,  1     an  values  of 
B  for  each  bed  . 
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and  B  had  to  be  accomplished.     The  first  step  was  to  equate 

P  u  2 
4£  »  f    '       °  (12) 

and,  from  equation  (9), 

f  -  lU± (1  ♦  0.6S4/Sd)2  +  B  ~  W) 

^u0  v  P 

Next,   a  Reynold's  number  representation  -which  could  be  used  as  an  abscissa 
•was  determined  by  inspection. 

j^.w^U^    b {lh) 

S*  A  (1  +  0.6St/Sp)2 

The  residual  term  was  considered  to  be  the  friction  factor, 

F  -  g  or  (15) 

Bw 

F  -  i-  +  1  (16) 

Re 

Re  and  F  were  solved  for  all  of  the  information  that  had  been  used  to 
determine  values  for  A  and  B»  They  were  also  determined  for  transition 
flow  data  that  were  not  used  for  determining  A  and  B,  "f"  was  determined 
by  equation  (12)  and  Re  by  use  of  equation  (lli)«  Values  for  A  and  B  were 
determined  by  equations  (10)  and  (11), 

Plate  VIII  compares  the  actual  values  for  F  and  Re  to  the  relationship 
suggested  by  equation  (16),  The  few  beds  tested  by  Brownell  and  co-workers 
(3)  represented  about  one  third  of  all  data  when  the  individual  test  runs  were 
counted.  Actually,  they  only  tested  five  beds,  or  about  7  percent  of  the 
number  of  beds  considered.  Four  out  of  five  of  their  runs  were  excluded  from 
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the  graph  in  order  that  unjustifiable  weight  -would  not  be  given  to  them. 
The  actual  values  of  F  and  Re  that  were  plotted  are  included  in  the 

appendix. 

DISCISSION 
Factors  of  General  Importance 

This  investigation  has  shown  that  pressure  drop  through  a  wide  variety 
of  packing  materials  can  be  accurately  predicted  by  reference  to  three 
easily  determined  properties  of  a  packed  bed.  These  properties  aret 
total  packing  surface,  S  ;  total  packing  perimeter,  C_;  and  free  volume  of 
the  packed  zone,  V.  Reference  to  these  properties  eliminates  the  necessity 
for  considering  highly  complex  methods  of  correlation  or  vague  terms  such 
as  "nomunal  particle  diameter"  or  "normal  voids"  • 

The  scope  of  packed  bed  variables  which  has  been  studied  is  summarized 
in  Tables  7  and  8.  The  accuracy  of  prediction  of  pressure  drop  is  also 
summarized  in  these  tables.  The  columns  headed  "average  deviation"  show 
how  well  the  pressure  drop  information  for  each  type  of  packing  material  is 
centered  upon  the  predicted  value  while  the  columns  headed  "root  mean 
square  deviation"  illustrate  the  average  error  involved  in  predicting  pressure 
drop. 

A  less  obvious  advantage  of  this  correlation  lies  in  the  fact  that 
coefficient  terms  do  not  vary  widely  for  different  types  of  packing 
materials.  The  coefficient  for  laminar  flow,  A,  varies  from  5»56  to  13.36 
or  2.lj  fold.  The  coefficient  for  turbulent  flow,  B,  varies  from  0.25  to 
0.079  or  3.2  fold.  These  ranges  of  variation  include  sparsely  packed  beds 
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to  dense  beds  and  large  ratios  of  column  diameter  to  particle  diameter  as 
well  as  ratios  of  column  diameter  to  particle  diameter  as  low  as  3:1*  The 
recent  correlation  by  Brownell  and  co-workers  (3)  involves  variation  of 
coefficient  terms  in  the  order  of  magnitude  of  10  fold  or  more  and  is 
incapable  of  predicting  pressure  drop  through  wire  packing.  Most  other 
correlations  cannot  be  compared  in  this  respect  because  they  were  not 
supposed  to  be  general  in  nature. 

The  greatest  deviations  of  pressure  drop  from  the  predicted  value 
were  encountered  in  the  case  of  information  published  by  Blake  (1),  Blake 
measured  pressure  drop  across  the  entire  packed  zone.  He  stated  in  his 
paper  that  the  packing  support,  where  units  such  as  Raschig  rings  usually 
assume  undesirable  orientation,  may  have  caused  the  overall  pressure  drop 
to  be  somewhat  higher  than  the  value  which  should  be  expected.  Estimations 
for  the  coefficient  for  turbulent  flow,  B,  scattered  most  widely  when  the 
ratio  of  column  diameter  to  particle  diameter  was  small.  This  is  a  logical 
consequence  since  very  few  packing  units  were  required  for  the  small  columns 
and  the  probability  of  any  one  such  bed  producing  a  representative  pressure 
drop  should  be  small. 

The  average  accuracy  of  prediction  of  pressure  drop,  including  question- 
able results  such  as  those  of  Blake,  was  found  to  be  -15  percent  for  laminar 
flow  and  -25  peicent  for  turbulent  flow. 

The  convergence  of  this  correlation  upon  those  of  other  persons  is  best 
illustrated  by  the  fact  that  results  of  other  investigators  forms  the  basic 
body  of  information  upon  which  the  present  conclusions  are  based.  The 
writer  performed  experimental  tests  which  were  primarily  designed  for  illus- 
trating the  effect  of  extremes  in  packed  bed  variables.  The  effect  of 
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"loose  pack"  and  "close  pack",  which  -was  studied  by  Oman  and  "Watson  (11), 
had  been  previously  treated  by  empirical  corrections.  This  investigation 
produced  good  correlation  for  such  extremes  without  reference  to  the  method 
of  packing  the  bed  save  that  randomness  should  be  maintained* 

The  correlation  for  B*  as  a  function  of  w/m  does  not  show  agreement 
with  Stoke1 s  law  for  freely  falling  objects  in  turbulent  flow.  Stoke1 s 
law  asserts  that  B  should  be  constant  where  the  objects  are  highly  dis- 

At 

persed.  The  approximation  that  has  been  used,  B  •  B(w/m),  implies  that 

m 

B  becomes  very  small  for  sparsely  packed  beds.  Pressure  drop  was  not 
measured  for  the  range  where  B*  is  predicted  to  be  very  small;  the  wire 
rings  constituted  the  limit  in  this  direction. 

Orifice  l^le;;.^  ,  An  oxperir.ent  VSl  conducted  to  determine  I  reason 
for  the  lack  of  convergence  upon  Stoke* s  law.  The  total  energy  loss  due 
to  flow  about  a  falling  object  was  compared  to  the  total  energy  loss 
through  an  orifice  in  a  pipe  line.  Freely  falling  objects,  as  treated  by 
Stoke1 s  law,  are  widely  dispersed  while  in  a  packed  column  the  objects  are 
encountered  frequently  by  the  fluid  stream.  Consequently,  loss  through 
widely  separated  orifices  was  compared  to  pressure  loss  through  a  series  of 
orifices  spaced  about  one  orifice  diameter  apart.  Nine  orifices  were  used. 
It  was  found  that  the  coefficient  for  pressure  loss  through  any  one  of  the 
widely  dispersed  orifices  was  considerably  higher  than  the  coefficient  for 
any  one  of  the  closely  spaced  orifices.  The  results  of  these  tests  were 
as  follows: 

Case  #1,  orifices  spaced  U*9h   orifice  diameters  apart. 
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Case  $2t  orifices  spaced  0.875  orifice  diameters  apart. 

B*  -  _££°6  -  0.11,56 

In  each  case,  velocity  was  based  on  the  orifice  area#  The  subscript,  o, 
refers  to  the  fact  that  an  orifice  was  considered.  The  experimental  pro- 
cedure and  apparatus  are  described  in  the  appendix. 

Case  //I  agrees  with  the  head  loss  coefficient  for  a  single  orifice. 
Case  #2  shows  that  only  3h*9   percent  as  much  energy  is  lost  v?hen  the  orifices 
are  spaced  to  compare  with  conditions  within  the  packed  bed.  This  analogy 
shows  that  Stoke* s  law  may  not  be  applicable  to  packed  beds.  It  implies 
that  coefficients  in  the  order  of  those  for  turbulent  flow  through  ducts  may 
be  approached  in  packed  beds.  Ducts  offer  much  less  resistance  to  flow 
for  a  given  total  surface  than  do  suspended  objects. 

Orientation  3 tear  the  Column  1  all.  The  coefficient  for  pressure  loss 
in  turbulent  flow,  B*,  suffers  large  depressions  as  the  column  diameter  is 
decreased.  This  depression  is  greater  than  that  suggested  by  the  reduction 
in  w/m  which  results  from  the  fact  that  small  columns  produce  less  dense 
beds  than  do  larger  columns.  Some  cross-sections  of  packed  beds  were 
exposed  to  determine  whether  the  arrangement  of  packing  units  was  such 
that  the  fluid  should  encounter  less  resistance  near  the  wall.  Photographs 
of  the  sections  that  were  damaged  least  during  preparation  are  sho-wn  in 
Plate  IX,  All  of  the  sections  are  shown  in  the  appendix.  These  sections 
were  prepared  by  settling  the  packing  into  a  thin  cement  slurry,  allowing 
the  cement  to  harden,  and  then  sawing  the  hardened  mass  into  crocs-sections 
at  intervals  of  about  one  inch.  Inspection  of  these  photographs  shows  that 
the  beds  represent  typical  degrees  of  packing  density  and  that  the  units  near 
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the  wall  are  arranged  similarly  to  the  units  in  the  interior  of  the  bed. 
Thus,  the  reduction  in  B  for  small  columns  cannot  be  attributed  to  the 
exLstance  of  larger  space  for  passage  near  the  •wall.  The  column  -wall  pro- 
bably tends  to  reduce  the  intensity  of  turbulence  -within  the  bed  so  that 
less  energy  is  lost  by  the  fluid  stream. 

How  "Well  Must  a  Packed  Bed  Be  Defined 

Brownell  and  Katz  (2)  felt  that  porosity  of  the  packed  bed  should  be 
determined  with  very  delicate  precision,  their  correlation  required  special 
knowledge  of  porosity.  Precision  certainly  does  not  detract  from  the 
validity  of  results,  but  it  is  often  difficult  to  measure  certain  proper- 
ties, such  as  porosity,  with  a  great  deal  of  accuracy.  Analysis  of  the 
proposed  equation  shows  just  how  errors  in  measuring,  or  predicting, 
packed  bed  variables  should  affect  the  accuracy  of  predicting  pressure 
drop.  Certain  errors  may  originate  from  definition  of  fluid  or  empty 
column  variables.  Lack  of  randomness  within  small  beds  might  also  contri- 
bute errors,  these  factors  are  not  included  in  the  following  discussion. 
Such  factors  as  the so  cannot  be  isolated  for  the  general  case,  however, 
the  fact  that  they  might  exist  is  sufficient  to  induce  necessary  precaution 
in  cases  where  they  may  become  predominant. 

Equation  (16)  summarizes  the  method  of  predicting  pressure  drop  for 
a  range  of  packed  bed  variables  that  includes  all  of  the  extremes  that 
might  be  encountered  in  its  application. 
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fors     1    <  w/n    <  35  (l6) 

0    <  w/Dt  <  0,3 
0    <  3t/Sp<0.3^ 

oo<vAt  <* 

The  range  of  applicability  may  be  larger,  but  only  the  rrnge  of  certainty 
is  stated.  The  range  of  certainty  encompasses  all  of  the  obseivstions  that 
have  been  cited. 

Equation  (17)  is  identical  to  equation  (16)  except  that  the  naperian 
base,  e,  has  been  substituted  for  the  base  10,  and  "w"  and  %n  have  been 
resolved  into  their  component  factors. 

Differentiation  of  equation  (17)  produces  equations  (18)  and  (19) ,  v.'hich 
are  of  direct  value  in  estimating  the  errors  which  might  arise  from 
inaccurate  estimation  of  the  different  variables.  Errors  in  estimating 
A  P/L .  for  highly  laminar  flow  are  sumarized  by  equation  (18),  Fluid 
and  empty  column  variables  are  considered  subject  to  no  error, 

di£ 

~||  -  0.0790  S§   /dS^  -  dV_  -  dC^V    dSj,  -     d£  |t  dSji 

iflp        V         CpyT      Sp       *  V    -**  Sp  ^p  (H-.6St/St>) 

when  flow  ia  laminar  (18) 
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For  highly  turbulent  flow,  error  is  sunmarized  by  equation  (19) • 


Lb 


=  -)..07*-£S£  „  *2fi\  +  3dS£  ^  2dV  ^  dC^  (19) 

DtVSp    CpJ    Sp    V    Cp 


when  flow  is  turbulent 
The  error  contributed  by  each  variable,  E,  is  sunmarized  in  Table  9» 


Table  9«  Contribution  of  incorrect  evaluation  of  packed  bed  variables 
to  error  in  predicting  pressure  drop. 


Variable  :     Laminar  flow  J    Turbulent  flow 


V 

E-  -(2  ♦  0.079-)^ 
¥                           m    V 

E     --2  dV 
▼               V 

Sp 

Es  -/2  ♦  0.1582 
1                  ra 

\ 

Es  -  (3  -  lu07£-)  2t 

L  1.2St 

Sp(l  ♦  .65^)^ 

P 

CP 

Ec  -  -0.0792  fSfc 

■    cp 

dCn 

E„  -  (14.07=-.  -  1)  _£ 

c            »t           cp 

Suppose  that  none  of  the  packed  bed  variables  is  to  be  allowed  to  con- 
tribute more  than  5  percent  to  the  error  in  predicting  pressure  drop.  Table 
9  shows  that  porosity,  or  V  should  be  known  within  2.5  percent  of  its  true 
value  when  flow  is  turbulent;  when  flow  is  laminar,  porosity  may  be  known 
within  2.5  percent  for  small  values  of  w/m,  or  within  1.57  percent  when 
w/m  reaches  the  upper  limit  of  15*   hen  flow  is  turbulent,  total  packing 
surface,  Sp,  should  be  known  within  1.67  percent  when  w/D^  is  small  and 
within  2.81  percent  when  w/D^  reaches  the  up<  er  limit  of  0.3;  laminar  flow 
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requires  that  packing  surface  be  known  within  2.5  percent  when  both  w/m 
and  St/Sp  are  small,  within  2.92  percent  when  S-j/Sp  reaches  its  upper 
limit  of  0.35*  within  l.lU  percent  when  w/m  reaches  its  upper  limit  of  15>, 
and  within  1.23  percent  when  w/m  and  S-^/Sp  both  reach  their  limiting 
values.     For  turbulent  flow,  the  total  packing  perimeter,   Cp,  should  be 
known  within  $  percent  for  small  values  of  w/D^  and  within  7.6  percent 
when  w/D^  reaches  its  upper  limit  of  0.3;  laminar  flow  requires  no 
knowledge  of  c    when  w/m  is  small  and  requires  accuracy  of  I4.22  percent  #ien 
w/m  reaches  its  upper  limit  of  15.     Transition  flow  requires  intermediate 
degrees  of  accuracy  for  these  variables. 

For  turbulent  flow,   the  sum  of  possible  errors  for  all  three  packed 
bed  variables  may  be  as  large  as  9»17  percent  for  small  values  of  w/D+ 
and  12,91  percent  for  the  limiting  w/D^  of  0.3  if  it  is  desired  to  maintain 
the  predicted  pressure  drop  within  15  percent  of  the  true  value.     Laminar 
flow  requires  that  this  sum  of  errors  be  within  5*0  percent  plus  any  large 
error  in  CL  vhen  w/m  is  near  zero  and  within  7.12  percent  when  S^/Sp  and 
w/ra  reach  their  respective  limits  of  0.35  and  1$. 

Sp  must  be  known  with  the  greatest  degree  of  accuracy,  V  requires  less 
accuracy,   and  C    may  be  less  well  defined  than  either  of  the  others  when 
each  of  the  variables  is  expected  to  contribute  the  same  degree  of  accuracy 
to  the  predicted  pressure  drop.     In  general,  an  error  in  pressure  drop  of 
less  than  -8.5  percent  will  result  if  the  value  of  each  of  these  variables 
is  known  within  -1  percent. 
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Determination  of  Packing  Surface  From  Pressure  Drop 

Knowledge  of  the  surface  area  of  inegular  objects  is  often  desired. 
Certain  relations  between  surface  area  and  particle  size,  catalytic  activity, 
mass  transfer  rates,  or  ionic  activity  exist.  Often,  a  quantitative  measure 
of  surface  area  can  be  used  to  determine  when  a  pulverizing  operation  is 
satisfactorily  completed,  or  when  fibers  are  of  a  desired  fineness  or  tex- 
ture. 

Most  tests  for  surface  area  have  been  conducted  by  allowing  flow  through 
the  material  in  question  to  be  laminar.  Allowing  the  flow  to  be  transitional 
would  require  very  cumbersome  calculations,  and  truly  turbulent  flow  requires 
tremendous  pressure  drops  that  are  not  required  when  laminar  flow  is  main- 
tained. Thus,  it  is  valid  to  assume  that  any  surface  area  determination  will 
require  the  use  of  equation  (20)  and  not  the  complete  equation  for  pressure 
drop. 

AP/L.  a  &*  g^£(1^0.6St/S13)2(10)0'03U3w/m 

or    &,  »    O.Dvfef  pg    .  A(io)-0.017lWm      m  0#6s 

for  highly  laminar  flow  (20) 

V,  L^,  ^A  ,  UQ,  Ap,  g,  and  St  are  always  kno-rcn  when  surface  area  is 
sought.  A  typical  determination  of  V  is  from  the  density  of  the  material 
and  the  total  volume  of  the  sample.  S^  is  usually  negligible  but  can  be 
included  when  the  necessity  arises. 

w/m  is  the  only  term  which  must  be  approximated  when  equation  (20)  is 
used.  Fortunately,  w/m  can  be  predicted  quite  accurately  from  porosity 
when  only  the  general  nature  of  the  particles  of  the  material  in  question 
is  known.  Table  10  shows  how  w/m  can  be  approximated  for  many  common 


shapes  that  are  approximated  by  fibers,  dusts,  crystals,  or  sands » 
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Table  10*  w/m  from  porosity  for  some  coiwion  geometrical  shapes. 


t         w/m  as  a  function  of  porosity 


Object 


Sphere 

Cube 

Tetrahedron 

Octahedron 

Circular  fiber 

Square  fiber 

Triangular  fiber 

Circular  disc 

Square  plate 

Equilateral  triangular 
plate 

Ribbon 


6  (1  -  V/Vt)/(V/Vt) 

N 

6.788(1  -  V/Vt)/(V/Vt) 

6.36k(l  -  V/Vt)/(V/Vt) 

2«  (1  -  V/Vt)/(V/Vt) 

8(1  -  V/Vt)/(V/Vt) 

5.196(1  -  v/vt)/(v/vt) 

(Diara^thickness)   (1  -  V/Vt)/(V/V  ) 
n 

7l73(edge/thickness)(l  -  V/Vt)/(V/Vt) 
2(width/thickness)  (1  -  V/Vt)/(V/Vt) 


Uniformity  of  particle  size  or  degree  of  conformity  to  the  geometrical 
shapes  Hated  effect  very  little  the  estimated  value  for  w/m.     Ordinarily, 
a  mixture  of  sizes  will  possess  a  slightly  lower  value  for  w/m  than  rill 
the  parent  particle  of  uniform  size.     Only  0.039$  of  the  error  in  estimating 
w/m  manifests  itself  as  relative  error  in  calculating  the  final  surface  area, 
that  is,  a  discrepancy  of  -1.0  in  evaluating  w/m  results  in  -3.95  percent 
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uncertainty  in  estimating  surface  area.  Ordinarily,  errors  in  evaluating 
w/m  should  range  between  zero,  for  materials  of  known  shape  to  -0.5  for 
materials  of  irregular  shape. 

Hoffing  and  Lcckhart  (15)  presented  information  by  which  the  method 
of  determining  surface  area  might  be  tested.  The  surface  area  of  dlatoma- 
cious  earth  was  determined  by  both  Nitrogen  adsorption  and  permeability. 
They  used  Carman^  conclusions  to  determine  surface  area  by  permeability. 
These  results  are  as  follows} 


Information  reported 

Vol.  of  cake,  V^ 
Area  of  cylinder,  A-^ 
Inverse  flow  rate 
Pressure  drop,  AP 
Porosity,  V/Vt 
Viscosity,  /x 
Surface  area,  Sg/Vp 
(by  methods  of  Carman) 

Surface  area,  Sp/Vp 
by  nitrogen  adsorption 

Derived  information 

w/m,  assuming  circular 
fibers 


Air  permeability 

0,756  cm3 
0.378  cm3 
8,87  sec/ cm-* 
70I4  gm/cm2 
0.711 
0.00018^  poise 

7ii600  cm2/cm3 
78800  cia2/cm3 


2.51 


Yjgter  permeability 

5.90  cm;J 
1.77  cm2 
27.3  sec/cur 
70I4  gm/cm2 
0.8*49 
0.009h7  poise 

714300  cm2/cm3 
78800  cn^/cm3 


1.12 


Surface  area,  Sp/V  •       -  o       1 

by  equation  (20)  75800  cray  cm3  71200  anycm-* 

This  example  does  not  show  the  full  advantage  of  equation  (20)   over 

that  proposed  by  Carman,  although  slightly  better  agreement  was  obtained 

by  equation  (20).     Equation  (20)   can  be  expected  to  apply  to  plate-like  or 

ribbon-like  materials  and  Carman's  equation  is  known  to  be  inapplicable  for 

such  materials. 
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One  Way  in  Which  Active  Surface  Might  be  Estimated 

Studies  concerned  with  mass  transfer  and  catalysis  in  packed  beds 
have  shown  that  all  of  the  surface  area  of  solid  particles  is  not  exposed 
to  the  transient  fluid  mass. 

The  proportion  of  unavailable  surface  might  be  comparable  to  an 
apparent  unavailable  free  volume.     The  following  equation,   similar  to  that 
of  Carman,  represents  how  pressure  drop  TOuld  be  expressed  for  hypothetical 
packed  beds  where  the  fluid  path  is  not  obstructed* 

AP  .  Jg0   s"VoH$p2        (1  +  Q#6  s   ,c  )2 
Lt        9  V3  F 

Actually,  the  fluid  path  is  obstructed  and  the  proportionality  between  the 
above  equation  and  equation  (20)  may  be  a  measure  of  the  surface  about 
which  appreciable  quantities  of  material  flow.     This  residual  term  is  pro- 
posed to  have  the  following  interpretation: 

active  surface  ■     available  volume     .  (y  aa    /2Q\-0,03l(3w/m 
total  surface  free  volume 

Turbulent  flow  is  not  understood  well  enough  so  that  the  effect  of  increasing 
flow  rate  upon  available  surface  can  be  predicted.     Present  knowledge  does 
not  preclude  the  possibility  that  the  same  fractional  free  volume  is 
available  in  turbulent  flow. 

Comparison  of  transfer  rates  in  packed  beds  requires  that  the  Reynold* s 
number  be  the  t;ame  for  each  bed.     Reynold's  number  has  been  determined  as 
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Re  -  I    3L£_k2 

A  yu(i  ♦  0.6  :>t/rp)2 

¥^w/200At/x  (i  ♦  0.6  st/rp)2(io)1-766w/Dt  *  o.03l,Vm 

Taecker  and  Hougen  (16)  determined  heat  transfer  coefficients  for 
several  packed  beds.  The  material  from  which  the  packing  was  made  and  the 
operating  conditions  were  maintained  very  nearly  constant  so  that  a  direct 
comparison  of  heat  transfer  coefficients  to  available  surface  can  be  made. 
The  following  information  illustrates  best  transfer  factors,  j^,  for  their 
beds  at  Re  m  3.0.  Values  of  w,  m,  snd  fL,  that  were  used  to  ascertain  the 
Reynold* s  number  are  also  included.  The  coefficients  axe  not  point  values 
but  have  been  obtained  from  plots  of  j^  vs.  Re  for  the  ran^e  of  interest. 

Packing        r^VA         m       w       &     (Vt/V)  (lo)"0^^/" 

0.0057  0.0350  0.103  0.976 

0.0057  0.0350  0.091  0.976 

0.0123  0.0698  0.095  0.396 

0.02ii3  0.1396  0.082  0.895 

0.0186  O.UiiO  0.062  0.813 

0.0038*  0.0358  0.11  0.806 

Required  free  volume  based  on  manufacturers  informs- tion  for  similar 
pecking  as  presented  by  Perry  (12). 

I'jgure  16  illustrates  the  manner  in  which  j^  decreases  as  the  pro- 
posed measure  for  available  surface  decreases,     k  decrease  in 
(V^/V)   (io)""0«03li3w/n  represents  &  proposed  decrease  in  available  surface. 
This  comparison  suggests  that  the  surface  available  for  transfer  nay  be 
proportion.-:!  to  (Vt/V)    (lo)-°*°^3Vm# 


Raschig  ring 
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Fig.  16 ,  Heat  transfer  factor  versus  a  measure  of 

available  surface  area,  Fe  constant  at  3«0. 
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Efficient  Packing 

The  packed  bed  is  used  "where  large  surfaces  are  required  for  mass 
and  heat  transfer  operations.  Increased  pressure  drop  through  a  bed 
increases  the  cost  required  in  accomplishing  a  given  amount  of  transfer* 
Thus,  packing  efficiency  may  be  defined  as  follows: 

packing  ef f  ■»  T^e  of  trans,   income  per  unit  trans, 
power  cons.     cost  per  unit  power 

Heat  transfer  rates  parallel  mass  transfer  rates  so  that  efficiency 
may  be  evaluated  from  the  standpoint  of  heat  transfer. 

packing  eff  ,  ■  — £ (income  ratio) 

U0At  AP 

This  efficiency  is  a  dimensionless  quantity  if  consistent  energy  units 
are  used  throughout.  Accurate  evaluation  of  efficiency  requires  specific 
knowledge  of  heat  transfer  coefficients  and  operating  conditions  for  each 
operation.  However,  efficiency  may  be  indexed  by  reference  to  some  well 
known  relationships. 

Most  transfer  processes  are  accomplished  where  flow  is  turbulent. 

Certain  general  relationships  may  be  established  for  this  case.  They  are 

as  follows: 

The  heat  transfer  coefficient,  h,  may  be  estimated  by  the  results 
of  the  investigation  by  Taecker  and  Hougen  (16),  The  Reynold's 
number  that  they  used  has  been  adjusted  to  include  "wH  rather  than 
their  approximation  for  particle  diameter.  This  equation  repre- 
sents their  tests  with  Raschig  rings  and  partition  rings. 


h  -  0.7230^  u0(s4-°-667(£i!ow)  -W* 


?h 


Pressure  drop  in  turbulent  flow  may  be  estimated  by  the  results  of 
this  investigation. 

AP  -It/°U02(Vm)Agm 

for  highly  turbulent  flow  and  large  diameter  columns. 

Substituting  these  approximations  into  the  equation  for  packing 
efficiency  produces  the  following  relationship. 

Packing  eff .  -  (inc.  ratio)  2^*°'67C°'3*T       .(SpW 

'^0.26^,0.10.^  2 .11        ^l.il 

The  group  of  fluid  and  empty  column  variables,  income  ratio,  k,  C,  A  T, 
M   j  P  ,   U  ,  and  constant  factors,  are  controlled  independently  of  the  type 
of  packed  column  so  that  (Sp/V.j.)jir  -*-w  *^  becomes  a  measure  of  the  expected 
amount  of  transfer  per  unit  of  pumping  power  input.  Including  "wtt  to  the 
1.0  power  rather  than  the  l.ljl  power  will  not  alter  the  accuracy  of  this 
approximation  by  a  greater  degree  than  the  order  of  accuracy  involved  in 
substituting  MwM  for  the  particle  diameter  used  by  Taecker  and  Hougen.  Thus, 
this  equation  may  be  considered  an  index  to  the  usefulness  of  a  packing 
material  for  transfer  operations. 

Packing  index  -  m(Vvt)  -  (V/Vt) 

w/ra     (w/m) 

when  flow  is  turbulent 

Similar  analysis  of  heat  transfer  during  laminar  flow  produces  this  index 

to  the  packing  materials 

Packing  index  -  (V/Vt)  (10)~°*°3ij3  w/ia 
when  flow  is  laminar 
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These  developments  have  not  considered  the  availability  of  surface  area 
■which  has  been  previously  discussed.  In  a  more  detailed  analysis,  the 
indexes  derived  above  should  be  multiplied  by  the  fraction  of  available 
surface.  Inclusion  of  the  proposed  measure  for  available  surface  -would  pro- 
duce these  inaexesi 

(10) -0.031*3  w/m 

Packing  ind.  ■  — - for  turbulent  flow 

(w/ra) 

■  (10)  -0.0686  w/m  for  laminar  flow 

The  measure  of  available  surface  does  not  change  the  overall  picture  as 
to  the  effect  of  w/m  on  packing  efficiency. 

Generally  speaking,  it  is  desirable  to  use  packing  materials  that 
produce  low  values  for  w/m.  The  range  of  w/m  noted  for  several  packing 
materials  that  were  studied  are  listed  below.  The  more  desirable  units 
are  placed  at  the  top  of  the  list. 

Packing  material  Range  of  w/m 

Wire  ring  1.21-  1.39 
Berl  saddle  3.95-  6.33 

Raschig  ring  ii.28-  8.23 
Cylinder  6.37-  9.65 

Sphere  6.80-10.53 

Mixing  wire  rings  with  other  packing  materials  produces  favorable 
values  of  w/m. 
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CONCLUSIONS 

Pressure  drop  in  packed  beds  has  been  found  to  depend  on  three  proper- 
ties of  the  packing.  These  properties  are  the  void  space  -within  the  bed, 
the  total  surface  of  the  packing,  and  the  perimeter  of  the  packing.  The 
perimeter  represents  boundries  which  must  be  circumvented  by  the  fluid  in 
passing  through  the  bed  and  has  been  defined  as  the  locus  of  tangent  points 
to  the  packing  that  would  be  intersected  by  a  line  that  moved  throughout  the 
bed  remaining  oriented  parallel  to  the  column  wall. 

Other  factors  incluencing  pressure  drop  are  the  size  of  the  column, 
fluid  density,  fluid  viscosity,  and  the  superficial  velocity  of  the  fluid. 
Normal  degrees  of  roughness  of  the  packing  material  do  not  noticeably 
influence  pressure  drop.  For  laminar  flow,  three  fifths  of  the  column  wall 
tends  to  reduce  pressure  drop  in  turbulent  flow,  this  effect  is  determined 
by  the  ratio  of  packing  width  to  column  diameter.  Packing  width  is  defined 
as  total  packing  surface  divided  by  total  packing  perimeter. 

Pressure  drop  through  beds  of  widely  varying  properties  and  for  flow 
ranging  from  completely  laminar  to  highly  turbulent  can  be  expressed  math- 
ematically or  graphically  by  reference  to  the  variables  which  are  mentioned 
above.  Graphical  correlation  of  pressure  drop  requires  consideration  of  the 
following  factors  . 

1.  F   ■  /^(Re),  primary  representation. 

2.  Re  -  n/°  V0B/^l   A(l  ♦  0.6  Sj/Sp) ? 

3.  F       -  f/B(w/m)  / 

2 

lu     f       -   APgm/Lt/°U0 


Table  of  nomenclature  is  on  page  7g. 
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5.  B   ■  /oXw/D^.),  auxiliary  representation. 

6.  A   ■  /(w/m),  auxiliary  representation. 
Mathematical  correlation  requires  use  of  the  following  equations  i 

ZPAt  -  V  Uo(1  *  0.6st/Sp)2/g»2  +  B*/3  U02/gm 

a  -^(io)°-03h30w/m      *i5* 

9 
B*  -  O^SCw/mKlO)-1-766  V»t    +25* 


These  conclusions  have  been  derived  from  information  representing  eight 
different  packing  materials  and  a  mixture  comprised  of  four  of  them.     They 
have  also  been  proven  applicable  to  materials  such  as  diatomacious  earth. 
The  range  of  variables  that  was  studied  is  as  follows  x 

0.3hk      4  V/Vt      ^  93.5 

0.020       ^  S^/Sp     ^       0.305 

1.21        ^  w/m       4:  11.11 

0.017      ^  w/Dt      ^        0.282 

0.00008  ^  Re        ^  97.0 

A  Reynold's  number  equal  to  unity  indicates  the  center  of  the  transi- 
tion region  as  flow  varies  from  laminar  to  turbulent. 

The  term,  w/m,  has  been  identified  as  an  index  to  the  power  required 
by  different  packing  materials  when  a  given  rate  of  mass  or  heat  transfer 
is  desired.     Small  values  of  w/m  indicate  the  minimum  in  pumping  costs. 


-it- 
Table  of  nomenclature  is  on  page  78  • 
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TABLE  OF  NOMENCLATURE 


ratio  between  the  effective  flow  area  and  A^. 

laminar  flow  constant  -  AP^/L^y*  UQ(1  *  0.6  S^/Sp)  -  B  N, 

the  predicted  value  for  A. 


a 
A 

Am.l. 

A« 

B 

Bra.l. 

B* 

B« 

Bo* 

c 


d 
D 


8 
DV 

e 
E 


cross-sectional  area  of  the  column, 
A(l  ♦  0.6  S^Sp)2. 

B*/(V*)  . 

predicted  value  for  B# 

turbulent  flow  constant  ■    6PgVl*t/°  U  2  -  A»/N. 

B*/(l  ♦  0.6  St/Sp)2. 

constant  for  an  orifice  •  APg/  P  U  • 

heat  capacity,  or  perimeter  contributed  by  a  specific  packing  unit. 

total  perimeter  of  the  packing, 

temperature,  degrees  Centigrade. 

centimeter. 

differential  element. 

equivalent  hydraulic  diameter,  or  nominal  particle  diameter  for  a 
specific  packing  unit. 

nominal  diameter  of  a  packing  unit. 

diameter  of  a  sphere  having  the  same  volume  as  the  packing  unit. 

diameter  of  the  column. 

effective  nominal  diameter  of  a  packing  unit  corresponding  to  a 
given  void  fraction. 

height  of  an  element  of  surface  roughness,  or,  2.71828. 

relative  error. 

error  in  predicted  pressure  drop  due  to  an  error  in  perimeter 
de  t erminat  io  n . 


79 


f 

*d 

fl 
ft 

F 

°F 
g 

gm 

h 

AH 

k 
K 

lb 
L 

m 

mm 

n 

N 

AP 

r 


erroi*  in  predicted  pressure  drop  due  to  an  error  in  surface 
determination. 

error  in  predicted  pressure  drop  due  to  an  error  in  porosity 
determination, 

APgm/Lt^u/. 

friction  factor  for  dense  packing  arrangement, 

frictiou  factor  for  loose  packing  arrangement, 

foot, 

friction  factor  •  f/B  • 

temperature,  degrees  Fahrenheit, 

acceleration  due  to  gravity,  taken  as  32 ,li>  ft,  per  sq,  sec,  in 
the  Manhattan  area, 

gram, 

heat  transfer  coefficient, 

vertical  displacement, 

heat  transfer  factor 

arbitrary  constant,  or  thermal  conductivity, 

arbitrary  constant, 

pound, 

true  length  of  flow  path,  or  length  of  a  specific  packing  unit, 

length  of  the  packed  zone  for  which  AP  is  measured, 

estimated  hydraulic  radius  ■  V/S   , 

pressure,  millimeters  of  Mercury, 

number  of  items, 

m^°Uo6u  • 

pressure  loss  due  to  frictional  resistance  along  L^, 

ratio  between  true  length  of  flow  path  and  L^, 


2/3 


hAt  /C/AW-* 
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He 

CR 

*P 

sq. 

sec. 

st 
% 

V 


w 
z 

e 


Reynold's  number  -  m  P  V0B*//JL  A(l  +0,6  St/S  )    . 

temperature,   degrees  Rankine. 

surface  of  the  packing  unit. 

square  » 

surface  of  a  sphere  having  the  same  volume  as  the  packing  unit. 

second* 

total  surface  of  the  packing. 

surface  of  the  column  wall. 

thickness  of  a  specific  packing  unit. 

true  fluid  velocity. 

velocity  based  on  the  empty  column. 

volume  of  a  specific  packing  unit, 

volume  of  the  packing  unit. 

free  or  void  volume  within  the  packed  zone. 

total  volume  of  the  packing. 

V  ♦  V  ■  volume  of  the  empty  column. 

Sp/Cp  ■»  width  of  the  barrier  to  be  circumvented  by  the  fluid. 

mass  rate  of  flow. 

ratio  between  effective  hydraulic  radius  and  m. 

arbitrary  function. 

constant  arising  because  of  the  geometrical  nature  of  a  duct. 

special  case  for  q   . 

incremental  element. 

angle  of  orientation  of  an  element  of  surface. 

absolute  viscosity. 


81 


IT      « 

•  3.1U159 

P      « 

»  mass  density  of  the  fluid* 

*       ' 

■  function  of. 

<   - 

■  less  than* 

^   ' 

•  less  than  or  equal  to* 

f^ 

•  approximated  by* 

cC 

»  proportional  to. 

T"  ' 

■  square  root. 
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APPENDIX 


EXPLANATION  OF  PLATE  X 

Cross-sections  of  packed  beds. 

Packing  materials  shomi? 

0.792  inch  glass  ball. 
1.028  inch  Berl  saddle. 
1.032  inch  Raschig  ring. 
0.512  inch  Berl  saddle. 
0.522  inch  Raschig  ring. 
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PLATE     X 
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1.      llj..  1     ;cntar\;    f.'o  :criptio:'is   of    packing;  uiaterials 


Hef.         PaOxCir.j "oscription 

(1)        Glass   rii  diameter  «   O.269   inch;   length  - 

0.261   inch;    thickness   *  0.032   inch. 

diameter  ■   0.233   inch;   length  = 
0.219   inch;    thickness   ■  0.032   inch. 

diameter  -  0.394-  inch;   length  a 
0.375   inch;    thickness   ■  0.032   inch. 

diameter  »  O.404  inch;  length  ■ 
O.449  inch;  thickness  «  0.030  inch. 

Glay  raschig  ring     diameter  =  1.00  inch;  length  = 

1.03  inch;  thickness  =  0.L1].  inch. 

(3)  oth   nickel    saddle      diameter   a  0.1316   inch;    surface   ■ 

6.3S  *  10-^  sq  ft;  volume  ■  ?.( 


x  10-0  cu  ft. 

010c  foot;  surfac 

,61 


Rough  nickel  packing  diameter  ».  0.010c  foot;  surface  a 

6. 3 J  x  10-^  sq  ft;  volume  =  7  •84 
x  10'°  cu  ft. 


dlay  Deri  saddle      diameter  =  1.00  inch;  surface  ■ 

O.O34.3  sq  ft;  volume  a  1.35  x  10"4 
cu  ft. 

Jiay  Raschig  ring     diameter  a  1.00  inch;  length  a, 

1.0C  inch;  volume  ■  I.93  x  10~4 
cu  ft. 

{$)  onal  pri         diameter  a  0.1f7  en;  lengt*  =  O.4.G 

cm. 


(11)        Celite    cylinder  diameter  ■  0.26/    inch;    length 

0.3.4r  inch. 


= 


:.q  diameter*  O.3C5  inch;  length  3 

0.397  inch;  *   O.O836  in. 

Olay  Eerl  saddle      diameter  a  0.5  inch;  surface  * 

378  sq  ft  per  cu  ft  of  packing. 
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Re 
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Re 
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).(  01 

-70 

O.00034 
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eoi 

.1)0122 

2 

60 

33 

2220 

-.00044 

3.ptT 
{.•hi 

7.61 

k? 

i. 

. 

41 

1.63 

5.79 

40 

1.47 

39 

1.52 

3. 

;  ■-• 

1.47 

0.60 

37 

1.43 

. 

. 

26 

,32 

..00143 

1230 

0.00086 

2. 

1. 

3. 

76 

1. 

3.61 

. 
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■ 
• 

74 

1. 

7.29 

1.24 

10.1 

72 

1.19 

12. 

71 

1. 

13.5 

1.32 

.20 

• 

l  • 

. 

1. 

.33 

59 

0.  . 

.6 

>8 

. 

. 

0.96 

.6 

131 

319 

O.OJ206 

132 

.  0147 

133 

,00 

0.00074 

115 

1. 

3.13 

114 

1.23 

^.19 

113 

0  c 

J- . 

p. 61 

112 

1. 

.50 

in 

1.07 

9. 

1. 

.70 

109 

0.91 

.15 

0. 

12.6 

107 

.62 

.5 

there 

G 

Tw  U 

0.00222 

H 

c. 00175 

T 

11 

0.00094 

7 

.  2 

':.00330 

00202 

1230 

c. 00074 

?? 

126 

0.00030 

331 

.  3  090 

1. 

. 00135 

1. 

4.00 
.93 

1. 

4.20 
C.75 

?.2Q 

1.59 

02 

1.94 

61 

1.71 

12.' 

..0 

1. 

.5 

c. 

17.2 

1. 

7.' 

17 

.  ?6 

.10 

16 

11.5 

- 

■  • 

. 

14 

.  044. 

20.3 

13 

. 

P     Z 

12 

3. 

.' 

10 

0.81 

31.3 

11 

.79 

33.  6 

12c 

157 

.   '10 

129 

306 

0.00321 

130 

0.00169 

i; 

2.10 

. 

120 

1. 

.. 

119 

<.42 

(  .36 

1: 

1. 

11. 

117 

i;6l 

.3 

116 

•> 
1.5c 

.1 

10 

1. 

5. 

105 

1.71 

1 .52 

12.7 

103 

1. 

.0 

102 

1.47 

.0 

l. 

37.2 

100 

1. 

45. ; 

161 

c. 

3. 

160 

0.00 

6.36 

159 

. 

10. 

0 .  o5 

6.79 

135 

n 

10. 

131+ 

0. 

15.; 

rl      D 

0.00403 

t 

0. 

I, 
4- 

022 

. .00035 

saddle 

216 

,  ■  .05 

9.56 

51 

2.40 

, .  70 

5o 

/•45 
16.1 

2.13 

4£ 

.03 

12.7 

47 

2.1 

•41 

19.5 

2.37 

23.7 

44 

2.52 

25 .( 

43 

2.1 

. 

2^ 

1.31 

13.6 

24 

1.32 

. 

23 

1.29 

.. 

22 

1.2; 

32.2 

21 

l.< 

3 

20 

1.21 

•  > 

19 

.22 

56.0 

90 

199 

.     .79 

0.003  >1 

342 

0.00292 

70 
67 

.7 

.71 

12.1 

0. 

10.6 
14.2  • 

jj> 

0.G3 

.3 

.CO 

19. 

53 

.76 

23.5 
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Ref.  Packing  Run; 


Rg   Run;'  F 


Re 


Run  - 


Re 


A.ut .    berl 

52 

0.71 

26.6 

152 

0.73 

7.60 

151 

0.67 

■ 
10.1 

saddle 

1>'0 

0.66 

13.5 

14? 

O.67 

17. 

148 

0.64 

21.8 

147 

0.66 

15.0 

Iii6 

.60 

22.2 

14-5 

0.60 

30.0 

HA 

0.59 

36.1 

^3 

0.56 

,.5 

"        Raschlg  125 

280 

. 

409 

0.00339 

127 

567 

0.00244 

ring 

124 

1. 

7 . 

123 

1.72 
1.6l 

10. 1 

122 

1. 

14. 5 

> 

99 

1.59 

11.3 

96 

Ik. 
.93 

97 

l., 
1.16 

21.6 

i.i+5 

25. 

- 

15o 

1.21 

157 

9.77 

156 

1.26 

13.  C 

155 

1. 

.2 

1ft 

1.47 

21.7 

153 

1.51 

142 

1. 11 

14.1 

141 

1.10 

20.1 

ll|0 

1.13 

21+. 

139 

9.87 

33- 

136 

1.21 

42.0 

137 

166 

1.10 

53.3 

184 

330 

0.00402 

165 

246 

0.00503 

265 

0.00 

177 

1.50 

24.0 

174 

19c 

1.47 

25.7 

176 

1. 
267 

29.9 

175 

1.43 

30.3 

290 

0.00410 

197 

0.00453 
23.0 

196 

I.4.C 

3?r 

O.OO46O 

193 

1.44 

23.4 

195 

I.46 

194 

30.3, 

192 

1.50 

30.4 

I69 

404 

0.00250  190 

0 . 00260 

191 

2340 

0.00042 

162 

1.39 

5.77 
9-69 

173 

1.34 

5.93 

172 

1.26 

9.04 

163 

1.22 

171 

1.21 

I3.4 

164 
166 

1.17 

13.9 

165 

1.11 

19.0 

170 

1.16 

19. 

1.10 

27.  : 

169 

1.10 

27. 

160 

1.06 

34  O 

167 

1.09 

35.2 

Mixture 

187 

7S3 

o. 00149 

745 

0.00149 

I.46 

5.86 

162 

1.35 

10.5 

1Q 

1.32 

13.1 

160 

1.26 

16.5 

179 

1.21 

29.1 

170 

1.19 

43..  5 

(1)      Raachig 

—— 

2.50 

1.96 

-»• 

2.10 

3.22 

•«• 

1.90 

k>52 

rin;5 

— 

1.76 

6.13 

•  •• 

1.66 

8.40 

-- 

1.56 

13.2 

wmmm 

2.12 

2.30 

— 

1.84 

1.56 

3.48 

-- 

1.76 

4.63 

•- 

1.6k 

5.74 

-- 

.01 

*»  ••- 

1.50 

10.4 

-- 

1. 

14.0 

-- 

2.05 

?• 

-- 

1.2? 

2. 

-- 

i.44~ 

.20 

— 

1.36 

-- 

7.55 

-- 

1.2k 

9.63 

-- 

1.19 

12. 

-- 

1.16 

15.9 

-- 

2.12 

4. 08 

*•  »• 

1.99 

j   'j'' 
^./c 

— 

1. 

7.70 

_   «• 

1.75 
i.5k 

10.2 

-- 

1.71 

12.5 

-~ 

1.61 

.6 

-- 

6.62 

-- 

1.12 

12.1 

-- 

0. 

16.0 

— 

0.94 

29.5 

— 

0.92 

39.3 

-- 

0.92 

h7-4 

(if)      Sphere 

•»  » 

3i.5 

0.0245 

__ 

20.5 

0.0415 
0.163 

_<•> 

9.35 

0.103 

*»  * 

7.55 

0.125 

-- 

5. 

— 

5.15 

0.20 

-- 

4.4 

?3 

M 

3.65. 

0.275 

— 

3.15 

0.355 

—  *» 

2.7 

0.46 



2.$ 

.  35 

— 

2.35 

0.595 

-- 

2.1 

0. 

-- 

2.0 

0.775 

-- 

1.95 

• .  hs 

-- 

1.59 

1.21 

-- 

1.45 

19.. 

3,*? 

**<m 

1.33 

2.3 

— 

1.17 

2.66 

«*     •• 

0.046    . 
0.265 

-- 

13.6 

0.073 



5.05 

0.225 

— 

4.6- 



3.05 

0.39 

U     15  (co;i  . 
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.'.  Pac  :in:    Ru 


He 


Run? 


F 


Re 


Run  ''■      F 


Re 


(4)  Sphere 


(3) 


Smooth 
saddle 


Rou 

saddle 


163 
180 

27 

101 

116 

81 

92 


2.3 
1.12 

\\  .0 

1. 
22. > 

3.2 
21+ 

2.7 

1.51 

1.3  3 
1.25 

U.5 

3.3 
12. 

?.35 

1. 

1. 
1.58 

1.39 

1.82 

1.13 

3390 
461 

•4 
35.2 
.27 

2.72 
2.07 
1.03 
O.96 


0.6?  5 

2.91 
0.2 

l.Ou 

o.c 

0.375 

0 .  0L5 

0.665 
3.0 

.2 

0.091 

O.38 

0.0 

0. 
1.68 
6.15 
1.05 
3.0 

0.133 
O.69 

3.3 

0.00029 
0.00220 
0.0108 
-.    22^ 
0.191 

.  ;55 
1.13 
.05 

28. 6 


—  lJj.6 

—  11.0 

—  3.5 

—  1.34 
--  10.1 

—  2.0 
--   11.7 

—  2.25 

—  i.5i 

—  1.33 

—  1.29 
--  6.0 

—  2.3 

—  c.55 

—  2.35 
--    1.60 

—  4-25 

—  1.90 
--    1.30 

—  4-3 

—  1.4.7 

—  0.91 

153  1130 
16^   ij.00 
164  65.2 
33  12.4 

4.27 

2.3o 

1. 

0.90 


106 

121 

71 


c6 


1.70 

0.007 

0.31 

1.94 
0.102 

0.61 
0.095 

0.8b  . 
3 .4- 

6 . 

0.197 
0.69 
0.131 
1.16 

2.15 
0.32 

1.03 

3.95 
0.235 

I.46 

8.8 

0.00077 
0.00255 
0.0140 

0.0735 
0.271 

0.714 

2.25 

13.3 


1.33 
7.8 

3.0 

i.oe 

6.0 

1.78 

6.35 
i.55 

1.49 
1.29 

33.5 
3.5 
2.15 

6.35 
2.2 

2.03 
2.1 

16.-2 
2.0 
1.25 


158   708 
175   126 
189  ^^ 
36  8.00 


111 

126 

76 

91 


3.23 
2.11 
1.23 
0.92 


2.10 
0.120 
0.405 
4.2 
0.168 
0.92 
0.190 
2.4 
3.65 
5.8 
0.026 
0.35 
0.765 
0.1?8 

1.43 
.05 
O.c 

2.35 
0.0625 
0.64 

2.35 


0.00141 
0.00724 

0.0164 

0.122 
0.401 
O.915 

445 


26.! 


136  10100  0, 


126 


1. 0001k  I4I      1560  0. OOO63   146        965   0.0012 
5i+o  0.00246  156       265  0.00465  161  98.9     0.0116 


171  52.7     0.0202 
10.3     c.120 
3.70  0.417 
2.5o  1.00 
1.4I    •  4.43 
0.52    ?b.? 
0.06    97.0 


212 
59 

103 

295 

277  o.c 

249    2.51  0.431 

202      1.26        3.01 


176  2c 6     0.0397 
191     7.80  0.169 
205     3.00  0.664 
76     1.95        1.28 
91     1.14       7.60 

51    0.75.  4?4 


181  17.4    0.0685 
196    5>.2$  0.264 


73 

61 

96 
56 


2.64  0.-C39 

1.56   2.53 
1.01  12. 
0.68  62.0 


4690  0.00016  93  2300  0.00034  98  1380  0.00062 
1030  0.00C94  106   520  c. 00166  113   302  0.00276 
4  0.00692  303  5o.6  0.0171  306  24.4  0.0362 
*  6.10  0.151  266  3.37  0.:: 

197  1.88  0.869   199  1.54  1.56 
206  1.94    .21  213  0.65  14.0 
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Rei".    ?aakin,-_-  Rua,'       :■• Re  Run.#      F ^e  Run-:'      ^  Re 

(3)    F  21:       0.75     214       231     0.66     I4.I.O       236     O.60     76.7 

8 ad die 
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Berl           21 

. 

0.00109 

22 

312 

0.00212 

19 

236 

0.002  ' 

saddle        2I4. 

216 

0.0032s 

25 

174 

0.00I4.03 

7 

.    149 

4?- 

0.00520 

27 

1C1, 

72 

3fr 

72.5 

0.0104 

35 

0.0161 

5 

.€ 

0.021 

4 

25.0 

0.0310 

60 

1". 

.  -487 

70 

9-5o 

0.0 

71 

6.65 

0.109 

76 

^ 

0.162 

124 

p 

472 

1.90 

0.706 

117 

141 

1.34 

113 

: 

2. 

109 

1.07 

3.02 

93 

0.72 

4.00 

5.< 

99 

0.C3 

7.05 

153 

O.63 
0.7ft 

9. 

151 

0.L3 

12.  c 

l'r6 

0. 

16 .9 

127 

26.0 

132 

0.73 

37.o 

137 

0.72 

52. 

142 

C.67 

76  4 

Ras  :  . . 

3I4 

.00147 

157 

224 

0 .  00429 

159 

125 

0.00746 

rln            103 

544 

0.0160 

IO4 

23.6 

0.0356 

ie.2 

0.0491 

105 

12. . 

0.0690 
0.1 

I4O 

9.25 

0.0 

lkl 

82 

7.24 

0.132 

6. 

122 

5.66 

0.162 

..57 

0.232 

llq. 

'-t  •  *  - 

0.2 

117 

3.65 

O.344 

120 

\  43 

0.392 

55 

241 

3.655 

90 

2.10 

0.81*5. 

92 

1.97 

1.61 

1.01 

5° 

1.91 

1.22 

4° 

1.81 

1. 

5  a 

1.70 

60 
37 

1.60 

1. 

2. 

% 

1. 

1.. 

2.21 
2. 

63 
73 

1.26 

1.33 

2.4O 
4.50 

70 

1.2: 

. 

21 

1.12 

7.70 

31 

1.12 

.  4' 

7 

1.07 

13.1 

17 

1.12 

17.. 

12 

0. 

21.6 

4 

1.03 

31.5 

Cylinder     1 

1 .  09 

16 . 0 

2 

1. 

17.2 

3 

1.06 

19.1 

4 

1. 

•  7 

.6 

1.4 

12.7 

1.10 

13.7 

•-> 
1 

1. 

1.11 

12. 

9 

1.03 

124 

23.6 

1.10 

13.0 

1. 

13.  c 

12 

0.90 

13 

0. 

27.0 

14 

O.67 

31.. 

0. 

31.2 

1: 

. 

17 

0.90 

21.1 

18 

0.90 

2^.0 

0.90 

.3 

20 

0. 

19.1 

21 

C.94 

204 

. 

. 

Ras                 1 

1.12 

i.e. 

31.0 
.0 

2 

> 

1.09 
1.09 

.2 

3^.6 

3 

6 

1.09 
1.09 

38.5 

2  1.6 

■  j> 
1 

1.10 

3?. 

8 

1.06 

9 

1.4 

23.8 

10 

:  .13 

26.0 

11 

1.12 

2{  . 

12 

1.12 

30.0 

1.16 

22.0 

4 

1.15 

23. s 

15 

]  .14 

2> 

1.21 

13.9 

17 

1.24 

15.2 

1.22 

16. 

1. 

.3 

20 

l.l6 

194 

21 

1.17 

20.1 

22 

1.15 

21. 

2? 

1.: 

23. c 

24 

0.86 

36, 

21 

0. 

0. 

51.6 

27 

0.94 

3if. 

. 

4l. 

29 

0. 

4J.0 

30 

c. 

30.0 

31 

.6 

32 

0. 

36 . 

33 

0. 

lj.1.2 

0. 

o.93 

22. 

36 

0.92 

25.7 

.91 

3o 

0. 

32.0 

39 

0.88 

37.3 
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R< 


Re 


Re 


(11)  Berl 

I 

1.09 

.0 
.2 

2 

1. 

..0 

s 

1.06 

• 

sadcle 

4* 

1.0? 

1.11 

.5 

1. 

-9 

. 

1.1 

• 

i. 

.  : 

9 

1. 

Ur  '  • 

10 

I. 

• 

11 

1.12 

23  .4 

12 

1.11 

32.3 

13 

1.10 

. 

11, 

1.09 

43.0 

15 

1.024. 

1.0 
.9 

16 

1.02 

J.  3 

17 

1.02 

57. 0 

If 

1.02 

19 

1. 

•  7 
1.6 

20 

1.00 

534 

21 

1.01 

.  i 

22 

1. 

23 

1.02 

50.4 
35.6 

2hr 

1.01 

57- 

25 

1. 

. 

1.03 

21 

1.02 

1*1.5 

26 

1.02 

.5 

29 

1. 01 

^0.0 

30 

1.01 

•  ' 

Sphere 

1 

1.014. 

2 

0. 

4 

3 

0.96 

22.1 

4 

1.03 

.1 

5 

1.02 

15.  t 

w 

1.01 

16.;; 

'•J 

1 

1. 

12.2 

1.07 
1.16 

13.0 

9 

1.0if 

.3 

10 

1.17 

7.69 

11 

e.77 

12 

1. 

.75 

\l 

.12 

10.3 

2:J. 

14 

1.11 

11.2 

15 

1. 

i  >.0 

1.12 

17 

1.: 

»i 

1. 

20.1 

19 

1  • 

22.1 

20 

1. 

.7 

21 

1. 

16.6 

1. 07 

. 

23 

1. 

23.2 

24 

X  .  U  : 

llwb 

i.G.+ 

15. 

2u 

1.02 

.2 

27 

1. 

16.5 

26 

1. 

2? 

0.99 

19.2 

30 

0.99 

21.5 

calculated 


the 


predicted  :.      of  A  and  '  .      Thus  the  accuracy  for  pre- 
dicting pressure   .        ■     srtained  by  corr.Dar.in--  F  to 
(1  -r  l/Re). 

Information  lackii  *s  is  tr<      in  the  same 

order  as  it  1   '  -G  12  or  13 . 
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Pressure  crop  through  nine  regularly  spaced 
one-i         i  orifices  in  a  three-fourth 

inch  pipe. 


Run  »  '-'   •     •  A.  H,  Inches 

i  ;. 

2 

T! 

n 

4- 

c 

7 
6 


inches 

Flow,     ^m/sec 

bo 

20.50 

36.9 

0.1^17 
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:0  flow  system  is  described  in  Piute  XI. 
v-  eJ  «  3A'i/9^i  all  terras  in  consistent  units  and  U  based  on 
the  area  of  the  orifice. 


EXPLANATION  OF  PLATE  XI 


A.  General  layout,  one  eighth  inch  ■  one  inch, 

a.  1'iater  Level  in  the  reservoir, 

b#  Thirty  gallon  tank  open  to  the  atmosphere. 

c.  1"  outlet  from  tank. 

d.  Internal  adaption  to  1  l/kn   pipe. 

e.  1  1/ii"  gate  valve,  wide  open  for  all  runs. 

f .  1  l/U"  close  nipple. 

g.  1  1/24"  elbow. 

h.  1  l/li"  close  nipple. 

1.  1  ]/!;"  x  ln  reducing  coupling. 

j.  1"  pipe  28  1/2"  in  length  exiting  to  the  atmosphere. 

AH.  Elevation  representing  energy  lost  in  the  orifices  . 

B.  The  dispersed  arrangement  for  the  orifices,  one-half  inch  ■  one  inch. 

k.  Upstream  insert  of  3/li"  pipe,  13  11/16"  long. 

1.  Downstream  insert  of  3/h"  pipe,  10  1/2"  long, 

m.  Long  spacer  nado  of  3A"  pipe,  1  15/61jM  long, 

n.  Orifice  made  of  washer  1"  o.d.  x  3/1;"  i.d.  x  3/l6"  thick. 

C.  The  close  arrangement  for  the  orifices,  one-half  inch  ■  one  inch. 

o.     Short  spacer  made  of  3/h"  pipe,   7/32"  long. 
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Pressure  drop  in  laminar  or  turbulent  flow  is  correlated  by  considering 
the  surface  of  the  packing  material,  the  porosity  of  the  jacked  bed,  the 
fluid  variables,  the  size  of  the  confining  column,  and  the  degree  to  which 
different  packing  materials  obstruct  flow.  The  effect  of  each  of  these 
factors  is  determined  according  to  a  simple  mathematical  or  graphical  solu- 
tion. 

Experimental  results  of  seven  investigators,  including  the  author,  are 
used  to  support  the  conclusions.  These  results  represent  seventy-five  beds 
packed  with  eight  different  types  of  packing  materials  and  a  mixture  of  four 
of  them.  The  correlation  is  of  such  a  nature  that  it  may  be  extended  to 
novel  packing  units. 

The  observed  accuracy  of  predicting  pressure  drop  is  l£  percent  for  lami- 
nar flow  and  2£  percent  for  turbulent  flow.  Accuracy  is  intermediate  to  these 
figures  for  transition  flow. 

Determination  of  surface  area  of  porous  media  by  permeability  is  discussed. 
Calculations  for  surface  area  are  outlined.  The  maximum  error  to  be  encountered 
in  determining  surface  area  is  estimated  to  lie  within  zero  to  2  percent  in 
cases  where  accurate  measurements  of  porosity,  pressure  loss,  and  flow  rate 
may  be  made,  and  where  the  approximate  shape  of  the  granules  is  known. 

Auxiliary  investigations  were  concerned  with  the  appearance  of  the  cross- 
section  of  a  packed  bed  and  the  effect  of  spacing  upon  pressure  loss  through 
a  series  of  orifices. 


